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We give a detailed account of the one-way quantum computer, a scheme of quantum computation 
that consists entirely of one-qubit measurements on a particular class of entangled states, the cluster 
states. We prove its universality, describe why its underlying computational model is different from 
the network model of quantum computation and relate quantum algorithms to mathematical graphs. 
Further we investigate the scaling of required resources and give a number of examples for circuits 
of practical interest such as the circuit for quantum Fourier transformation and for the quantum 
adder. Finally, we describe computation with clusters of finite size. 
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I. INTRODUCTION 

Recently, we introduced the scheme of the one-way 
quantum computer (QC C ) 0. This scheme uses a given 
entangled state, the so-called cluster state Q, as its cen- 
tral physical resource. The entire quantum computation 
consists only of a sequence of one-qubit projective mea- 
surements on this entangled state. Thus it uses measure- 
ments as the central tool to drive a computation - @ • 
We called this scheme the "one-way quantum computer" 
since the entanglement in the cluster state is destroyed 
by the one-qubit measurements and therefore it can only 
be used once. To emphasize the importance of the clus- 
ter state for the scheme, we use the abbreviation QC C for 
"one-way quantum computer" . 

The QC C is universal since any unitary quantum logic 
network can be simulated on it efficiently. The QC C can 
thus be explained as a simulator of quantum logic net- 
works. However, the computational model that emerges 
for the QC C makes no reference to the concept of 
unitary evolution and it shall be pointed out from the 
beginning that the network model does not provide the 
most suitable description for the QC C . Nevertheless, the 
network model is the most widely used form of describing 
a quantum computer and therefore the relation between 
the network model and the QC C must be clarified. 

The purpose of this paper is threefold. First, it is to 
give the proof for universality of the QC C ; second, to 
relate quantum algorithms to graphs; and third, to pro- 
vide a number of examples for QC c -circuits which are 
characteristic and of practical interest. 

In Section [n] we give the universality proof for the de- 
scribed scheme of computation in a complete and detailed 
form. The proof has already been presented to a large 
part in Q. What was not contained in was the ex- 
planation of why and how the gate simulations on the 
QC C work. This omission seemed in order since the im- 
plementation of the gates discussed there (CNOT and 
arbitrary rotations) require only small clusters such that 
the functioning of the gates can be easily verified in a 
computer simulation. For the examples of gates and sub- 
circuits given in Section IIVI this is no longer the case. 
Generally, we want an analytic explanation for the func- 
tioning of the gate simulations on the QC C . This expla- 



nation is given in Section |HFj and applied to the gates of 
a universal set in Section III Gl as well as to more compli- 
cated examples in Section llVl 

In Section III HI we discuss the spatial, temporal and 
operational resources required in QC c -computations in 
relation to the resources needed for the corresponding 
quantum logic networks. We find that overheads are 
at most polynomial. But there do not always need to 
be overheads. For example, as shown in Section III II 
all QC c -circuits in the Clifford group have unit logical 
depth. 

In Section IIIII we discuss non-network aspects of the 
QC C . In Section Hit Al we state the reasons why the net- 
work model is not adequate to describe the QC C in every 
respect. The network model is abandoned and replaced 
by a more appropriate model • This model is described 
very briefly. 

In Section IIII Bl we relate algorithms to graphs. We 
show that from every algorithm its Clifford part can be 
removed. The required algorithm-specific non-universal 
quantum resource to run the remainder of the quantum 
algorithm on the QC C is then a graph state 0- All that 
remains of the Clifford part is a mathematical graph spec- 
ifying this graph state. 

In Section lTVl we give examples of larger gates and sub- 
circuits which may be of practical relevance, among them 
the QC c -circuit for quantum Fourier transformation and 
for the 77-qubit adder. 

In Section we discuss the QC C computations on fi- 
nite (small) clusters and in the presence of decoherence. 
We describe a variant of the scheme consisting of re- 
peated steps of (re-)entangling a cluster via the Ising 
interaction, alternating with rounds of one-qubit mea- 
surements. Using this modified scheme it is possible to 
split long computations such that they fit piecewise on a 
small cluster. 



II. UNIVERSALITY OF QUANTUM 
COMPUTATION VIA 
ONE-QUBIT-MEASUREMENTS 



In this section we prove that the QC C is a universal 
quantum computer. The technique to accomplish this is 
to show that any quantum logic network can be simulated 
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efficiently on the QC C . Before we go into the details, let 
us state the general picture. 

For the one-way quantum computer, the entire re- 
source for the quantum computation is provided initially 
in the form of a specific entangled state -the cluster state 
0- of a large number of qubits. Information is then 
written onto the cluster, processed, and read out from 
the cluster by one-particle measurements only. The en- 
tangled state of the cluster thereby serves as a universal 
"substrate" for any quantum computation. It provides in 
advance all entanglement that is involved in the subse- 
quent quantum computation. Cluster states can be cre- 
ated efficiently in any system with a quantum Ising-type 
interaction (at very low temperatures) between two-state 
particles in a lattice configuration. 

It is important to realize here that information pro- 
cessing is possible even though the result of every mea- 
surement in any direction of the Bloch sphere is com- 
pletely random. The mathematical expression for the 
randomness of the measurement results is that the re- 
duced density operator for each qubit in the cluster state 
is hi. The individual measurement results arc random 
but correlated, and these correlations enable quantum 
computation on the QC C . 

For clarity, let us emphasize that in the scheme of the 
QC C we distinguish between cluster qubits on C which 
are measured in the process of computation, and the log- 
ical qubits. The logical qubits constitute the quantum 
information being processed while the cluster qubits in 
the initial cluster state form an entanglement resource. 
Measurements of their individual one-qubit state drive 
the computation. 

information flow 
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FIG. 1: Simulation of a quantum logic network by measuring 
two-state particles on a lattice. Before the measurements the 
qubits are in the cluster state \(j>)c of Q. Circles symbolize 
measurements of a z , vertical arrows are measurements of er B) 
while tilted arrows refer to measurements in the x-y-plane. 



To process quantum information with this cluster, it 
suffices to measure its particles in a certain order and in a 
certain basis, as depicted in Fig.^ Quantum information 
is thereby propagated through the cluster and processed. 
Measurements of o^-observables effectively remove the 
respective lattice qubit from the cluster. Measurements 
in the a x - (and a v -) eigenbasis are used for "wires" , i.e. 



to propagate logical quantum bits through the cluster, 
and for the CNOT-gate between two logical qubits. Ob- 
servables of the form cos(y>) <j x ± sin(i^) a y are measured 
to realize arbitrary rotations of logical qubits. For these 
cluster qubits, the basis in which each of them is mea- 
sured depends on the results of preceding measurements. 
This introduces a temporal order in which the measure- 
ments have to be performed. The processing is finished 
once all qubits except a last one on each wire have been 
measured. The remaining unmeasured qubits form the 
quantum register which is now ready to be read out. At 
this point, the results of previous measurements deter- 
mine in which basis these "output" qubits need to be 
measured for the final readout, or if the readout mea- 
surements are in the a x -, cr y - or cr z -eigenbasis, how the 
readout measurements have to be interpreted. Without 
loss of generality, we assume in this paper that the read- 
out measurements are performed in the aveigenbasis. 



A. Cluster states and their quantum correlations 

Cluster states are pure quantum states of two-level sys- 
tems (qubits) located on a cluster C. This cluster is a 
connected subset of a simple cubic lattice Z d in d > 1 
dimensions. The cluster states |</ , {«}) c obey the set of 
eigenvalue equations 



^ (o) l^w) c = (-i) Kt 'l^w) C! 

with the correlation operators 



4 a) <g> 

fcGnbgh(a) 



(6) 



(1) 



(2) 



Therein, {n a <G {0, 1}| a G C} is a set of binary param- 
eters which specify the cluster state and nbgh(a) is the 
set of all neighboring lattice sites of a. All states |0{ K }) C 
are equally good for computation. A cluster state is com- 
pletely specified by the eigenvalue equations (I}. To see 
this, first note that two states |0{ K })c and |0{s})c which 
obey a set of equations JTJ but differ in at least one 
eigenvalue are orthogonal. This holds because if there 
exists an a € C such that, say, K^\4>{ K })c = \4>{k,})c 
and K^\cj) {k} )c = -|<£{«})c, then c(</>{s}|<£{ re })c = 

c (</>{«} |^ (Q) |0 {K })c - -c<0{*}|^{„}>c = 0. From the set 
of states which obey QJ with the eigenvalues specified 
by {k} a representative \4>{ K })c is taken. There are 2l c l 
such classes of states, and hence 2' c l mutually orthog- 
onal representatives \4>{ K })c- Therefore, the representa- 
tive cluster states form a basis {|</>{ K })c | {k} 6 {0, 1}' C '} 
of the |C|-qubit Hubert space. To that end, let us 

that obeys with the 



now consider a state 



same {k}o as \4>{ K } )c, an d expand it into the above ba- 
sis. One finds \<f>\ K}o )c = £{«} c<0{„}|^ {b}o }c |0{k})c = 
c{4>{ K } \<t>' {K}o )c I0{«} o )c- Hence, two states W {k}o )c and 
\4>{k} )c which obey with the same set {k} are the 
same modulo a possible global phase. Consequently, any 



3 



method that creates a state obeying equations with 
a specific set {n a \a <G C} creates the same state. 

The eigenvalue equations JU and the quantum corre- 
lations they imply are central for the described scheme of 
computation. Also, they represent a very compact way 
of characterizing the cluster states. To reflect this in the 
presentation, the discussion in this paper will be based 
entirely on these eigenvalue equations and we will never 
need to work out some cluster state in any specific ba- 
sis. In fact, to write down a cluster state in its explicit 
form would be quite space-consuming since the minimum 
number of required terms scales exponentially with the 
number of qubits Q , and for computation we will be go- 
ing to consider rather large cluster states. Nevertheless, 
for illustration we give a few examples of cluster states 
with a small number of qubits. The cluster states on a 
chain of 2, 3 and 4 qubits, fulfilling the eigenvalue equa- 
tions QJ with all K a = 0, are 

l^)c 2 = ^=(|0)i|+> 2 + |l)i|-) 2 ), 

\<t>)c 3 = ^|(l+>l|0> 2 |+>3 + |-)l|l>2|->3), 
\4>) Ci = ^l+>l|0>2|+> 3 |0)4 + i|+)l|0) 2 h) 3 |l>4, 

+ ^|-)l|l)2h)3|0)4 + i|-)l|l) 2 |+)3|l)4, 

(3) 

with the notations 

|0)„ := \0) a . z =a z a) \0) a . z , 
|1>« := |l)a,, = -ai a) |l) a ,„ (4) 
|±> a := ^(|0) o ±|l) a ). 

The state \4>}c 2 1S local unitary equivalent to a Bell state 
and \4>)e 3 to the Greenberger-Horne-Zeilinger (GHZ) 
state. \4>)ci is not equivalent to a 4-particle GHZ state. 
In particular, the entanglement in l</>)c 4 cannot be de- 
stroyed by a single local operation 

Ways to create a cluster state in principle arc to mea- 
sure all the correlation operators , a € C of on an 
arbitrary |C|-qubit state or to cool into the ground state 
of a Hamiltonian Hk = ~%Saec K aK^ a \ 

Another way -likely to be more suitable for realiza- 
tion in the lab- is as follows. First, a product state 
|+)c = ®aeC 1+)° ^ s prepared. Second, the unitary 
transformation S^ c ', 

5 (c) = Yl S ab , (5) 

a,b£C\b— a£7 d 

is applied to the state |+). Often we will write S 
in short for S^ C K In J5J, for the cases of dimension 
d = 1, 2, 3, we have 7l = {1}, 72 = {(1, 0) T , (0, 1) T } and 
73 = {(1,0,0) T ,(0, 1,0) T ,(0,0,1) T }, and the two-qubit 
transformation S ab is such that the state |1) Q ® |l)b ac- 
quires a phase of it under its action while the remaining 
states |0) a ® 1 0) t, , |0) Q <g> 1 1) t. and |1) Q ® |0)f, acquire no 
phase. Thus, S ab has the form 

S ab = \(l + o^+o^-a^®a^). (6) 



The state \+)c obviously obeys the eigenvalue equations 

o x \+)c = |+)c Va e C and thus the cluster state \<j)}c 
generated via S obeys 

\4>) c = So^Si \cj>)c , VaeC. (7) 
To obtain S<Jx S\ observe that 

S a b(T {«) S a b \ = .W^W 

and 

S ab a^S ab ^ =4 c) . VceC\{a,&}. (9) 

Further, the Pauli phase flip operators <j z commute with 
all S ab , i.e. 

S ab ai d) S ab] = ff M i VdeC. (10) 
Now, from ©, © and (JTUl) it follows that 

SvP&=<tM ® <7< 6 >. (11) 

6£nbgh(a) 

Thus, the state \(j))c generated from |+)c via the trans- 
formation S as defined in JHJ) does indeed obey eigenvalue 
equations of form , with 

K a = 0, MaeC. (12) 

Note that all operations S ab in S mutually commute 
and that they can therefore be carried out at the same 
time. Initial individual preparation of the cluster qubits 
in |+)aec can also be done in parallel. Thus, the creation 
of the cluster state is a two step process. The temporal 
resources to create the cluster state are constant in the 
size of the cluster. 

If a cluster state is created as described above this leads 
to the specific set of eigenvalues in Q specified by the 
parameters n a in (|12|) . As the eigenvalues are fixed in 
this case, we drop them in the notation for the cluster 
state \(j>)c ■ Cluster states specified by different sets {k q } 
can be obtained by applying Pauli phase flip operators 
<j z a \ To see this, note that 

a^K^a^ = (-l) s "<» JjfW. (13) 

Therefore, 

(g) |0 {Ka} ) c = |0 {Ka+ A Ka }>c (14) 

aec 

where the addition for the K a is modulo 2. 

The transformation S defined in JSJ is generated by 
the Hamiltonian 

1 rr (a) 1 rr {b) 

H = hg £ 1 £ , (15) 

a,b£C\b— a£f d 
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and S is of the form 



S = cxp — m 



l-ai a) l-aj b) 
2 2 



(16) 



Expanding the exponent in (|16f> . one obtains 



S 



n e ' 



exp 



7T 
Z — (7 

4 



(a) 



cxp 



V 4 z 



x cxp 



f-i E «w) 

y a,6GC|6-a67d / 



(17) 

We find that the interaction part Hi of the Hamiltonian 
H generating S is of Ising form, 



= E (18) 

a,b£C\b-a£y d 



and, since the local part -£q oca i of the Hamiltonian com- 
mutes with the Ising Hamiltonian Hi, the interaction S 
generated by H is local unitary equivalent to the unitary 
transformation generated by a Ising Hamiltonian. 

For matter of presentation, the interaction S ab in © 
and, correspondingly, the local part of the Hamiltonian H 
in (|15f) has been chosen in such a way that the eigenvalue 
equations JQ) take the particularly simple form with K a = 
for all a G C, irrespective of the shape of the cluster. 

Concerning the creation of states that are useful as a 
resource for the QC C , i.e. cluster- or local unitary equiv- 
alent states, all systems with a tunable Ising interaction 
and a local cr z -type Hamiltonian, i.e. with a Hamiltonian 



H> = J2*E a ^+h 9 -^ £ 



(19) 



a,b£C\b— aEya 



are suitable, provided the coupling g(t) can be switched 
between zero and at least one nonzero value. 

Even this condition can be relaxed. A permanent Ising 
interaction instead of a globally tunable one is sufficient, 
if the measurement process is much faster than the char- 
acteristic time scale for the Ising interaction, i.e. if the 
measurements are stroboscopic. If it takes the Ising in- 
teraction a time Xi s i n g to create a cluster state \<j))c 
from a product state |+)c, then the Ising interaction 
acting for a time 2Ti s i ng performs the identity operation, 
5'(C) J 5'(C) _ \<fi) - Therefore, starting with a product state 
\+)c at time t — evolving under permanent Ising inter- 
action, stroboscopic measurements may be performed at 
times (2k + 1) Ti sing , k G N. 

Some basic notions of graph theory will later, in the 
universality proof, simplify the formulation of our speci- 
fications. Therefore let us, at this point, establish a con- 
nection between quantum states such as the cluster state 
of JIJ and graphs. The treatment here follows that of 0, 
adapted to our notation. 



First let us recall the definition of a graph. A graph 
G(V, E) is a set V of vertices connected via edges e from 
the set E. The information of which vertex a G V is 
connected to which other vertex b G V is contained in 
a symmetric \V\ x \V\ matrix T, the adjacency matrix. 
The matrix T is such that T a b = 1 if two vertices a and b 
are connected via an edge e G E, and r a ;, = otherwise. 
We identify the cluster C with the vertices Vc of a graph, 
C = Vc, and in this way establish a connection to the 
notion introduced earlier. 

To relate graphs to quantum mechanics, the vertices 
of a graph can be identified with local quantum systems, 
in this case qubits, and the edges with two-particle in- 
teractions [HjGSj in the present case er z <7 z -interactions. 
If one initially prepares each individual qubit a in the 

fa) Ka 

state (<7 Z ) \+) a and subsequently switches on, for an 
appropriately chosen finite time span, the interaction 



Hg(v,b) = hg 51 

(a,b)eE 



1 - ai a) 1 



(20) 



with (a, b) G E denoting an edge between qubits a and 
b, then one obtains quantum states that are graph code 
words as introduced in 0. Henceforth we will refer to 
these graph code words as graph states and use them in a 
context different from coding. The graph states \4>{k})g 
arc defined by a set of eigenvalue equations which read 



4 a) (gj I - 

bev 



(6) 



\4>{K}) G = {-l) Ka \4>{K}) G , (21) 



with n a G {0, 1} Va G V. Here we use G instead of V as 
an index for the state \<j>) as the set E C V x V of edges 
is now independent and no longer implicitly specified by 
V as was the case in 

Note that cluster states Q arc a particular case of 
graph states (|21|l . The graph G(C,Ec) which describes 
a cluster state is that of a square lattice in 2D and that 
of a simple cubic lattice in 3D, i.e. the set Eq of edges is 
given by 



E c = {(a,6)|a,6GC,6Gnbgh(a)}. 



(22) 



Let us at the end of this section mention how clus- 
ter states may be created in practice. One possibility 
is via cold controlled collisions in optical lattices, as de- 
scribed in 0| . Cold atoms representing the qubits can 
be arranged on a two- or three dimensional lattice and 
state-dependent interaction phases may be acquired via 
cold collisions between neighboring atoms [l2l| or via tun- 
neling 01 Fot a suitable choice of the collision phases tp, 
(p = 7T mod 2-7T, the state resulting from a product state 
\+)c after interaction is a cluster state obeying the eigen- 
value equations lJT]l. with the set {n a , a G C} specified by 
the filling pattern of the lattice. 
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B. A universal set of quantum gates 

To provide something definite to discuss right from the 
beginning, we now give the procedures of how to realize 
a CNOT-gate and a general one-qubit rotation via one- 
qubit measurements on a cluster state. The explanation 
of why and how these gates work will be given in Sec- 
tion CEG1 

1 2 3 4 5 6 7 
control X 
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12 3 4 

H 
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1 2 3 4 5 
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1 2 3 4 5 

E)[x x □ 



z-rotation 
(e) 

1 2 3 4 5 

x x n x n 

7r/2-phase gate 



FIG. 2: Realization of elementary quantum gates on the QC C . 
Each square represents a lattice qubit. The squares in the 
extreme left column marked with white circles denote the in- 
put qubits, those in the right-most column denote the output 
qubits. 



A CNOT-gate can be realized on a cluster state of 15 
qubits, as shown in Fig. [21 All measurements can be 
performed simultaneously. The procedure to realize a 
CNOT-gate on a cluster with 15 qubits as displayed in 
Fig. His 

Procedure 1 Realization of a CNOT-gate acting on a 
two-qubit state IV'in)- 



1. Prepare the state 

|*m>C 15 = |^m)l,9 <8> (® ig e 15 \{l,9} 



2. Entangle the 15 qubits of the cluster C15 via the 
unitary operation 

3. Measure all qubits of C15 except for the output 
qubits 7, 15 (following the labeling in Fig. EJ - The 
measurements can be performed simultaneously. 
Qubits 1, 9, 10, 11, 13, 14 are measured in the o x - 
cigenbasis and qubits 2-6, 8, 12 in the er. y -cigcnbasis. 



Dependent on the measurement results, the following 
gate is thereby realized: 

U'cnot = U^cnot CNOT(c, t). (23) 
Therein the byproduct operator Uy,,cnot has the form 



CNOT 

(c) 
lx 

(*) 

71 C) 



(c)Tx C) (t)7< t) (c) 7< c) {t pi t} 
Ox Ox O z O z , with 



S2 + S3 + S 5 + S 6 

S2 + S3 + S 8 + Sio + Sl2 + s u 

Si + S3 + S 4 + S 5 + Sg + Sg + Su + 1 

sg + Su + S13. 

(24) 



Therein, the Sj represent the measurement outcomes s, 
on the qubits i. The expression l|24|l is modified if re- 
dundant cluster qubits are present and/or if the cluster 
state on which the CNOT gate is realized is specified by 
a set {« a } different from l|12f) . see Section lll CI This con- 
cludes the presentation of the CNOT gate, the proof of 
its functioning is given in Section lll CI 

An arbitrary rotation Unot G SU(2) can be realized 
on a chain of 5 qubits. Consider a rotation in its Euler 
representation 



U R ot^,V,C} = U x {C]U z [r,}U x ^}, 
where the rotations about the x- and z-axis are 



U x [a] 
UJa] 



exp 
exp 



(- ia f) 

( ■ a 

\ 2 



(25) 



(26) 



Initially, the first qubit is prepared in some state IV'in) > 
which is to be rotated, and the other qubits are prepared 
in |+). After the 5 qubits are entangled by the unitary 
transformation S, the state IV'in) can be rotated by mea- 
suring qubits 1 to 4. At the same time, the state is also 
swapped to site 5. The qubits 1 .. 4 are measured in ap- 
propriately chosen bases 



\0} j + e i ^\l) j lO^-e^ll), 



V2 



V2 



(27) 



whereby the measurement outcomes Sj € {0, 1} for j = 
1 ..4 are obtained. Here, Sj = means that qubit j is 
projected into the first state of Bj(ifj). In (|27|l the ba- 
sis states of all possible measurement bases lie on the 
equator of the Bloch sphere, i.e. on the intersection of 
the Bloch sphere with the x-y-plane. Therefore, the mea- 
surement basis for qubit j can be specified by a single pa- 
rameter, the measurement angle tpj. The measurement 
direction of qubit j is the vector on the Bloch sphere 
which corresponds to the first state in the measurement 
basis Bj(ifj). Thus, the measurement angle tpj is the an- 
gle between the measurement direction at qubit j and 
the positive x-axis. In summary, the procedure to realize 
an arbitrary rotation t/Rot[£, f], Q, specified by its Euler 
angles £, ?/, C, is this: 
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Procedure 2 Realization of general onc-qubit rotations 
U Rot £ SU(2). 

1. Prepare the state |* in )c 5 = |^in)i ® (®f =2 

2. Entangle the five qubits of the cluster C5 via the 
unitary operation 

3. Measure qubits 1 - 4 in the following order and basis 

measure qubit 1 in Si(0) 
02 measure qubit 2 in # 2 (-£(-l) Sl ) / 28 x 
I2J3 measure qubit 3 in £> 3 (—7/ (— l) S2 ) ^ ' 

04 measure qubit 4 in S 4 {-( (-l) Sl+S3 ) 

Via Procedure |21 the rotation U' Rot is realized: 

U' Rot [^v,C} = U SiRot U Rot [^r,,C}. (29) 
Therein, the random byproduct operator has the form 
C/E,i? ot = ^ 2+S4 < +S3 . (30) 

It can be corrected for at the end of the computation, as 
will be explained in Section Til El 

There is a subgroup of rotations for which the realiza- 
tion procedure is somewhat simpler than Procedure 
These rotations form the subgroup of local operations in 
the Clifford group. The Clifford group is the normalizcr 
of the Pauli group. 

Among these rotations are, for example, the Hadamard 
gate and the 7r/2-phase gate. These gates can be realized 
on a chain of 5 qubits in the following way: 

Procedure 3 Realization of a Hadamard- and tt/2- 
phase gate. 

1. Prepare the state |^ in )c 5 = |^in)i ® (<8>i=2 !+>*)• 

2. Entangle the five qubits of the cluster C5 via the 
unitary operation 

3. Measure qubits 1-4. This can be done simultane- 
ously. For the Hadamard gate, measure individu- 
ally the observables a x X \ &y , a< y\ a y '■ For the 

In 1 *. (!) ( 2 ) (3) (4) 

7r/2-phase gate measure a x , a x , a y , a x 

The difference with respect to Procedure |2 for general 
rotations is that in Procedure no measurement bases 
need to be adjusted according to previous measurement 
results and therefore the measurements can all be per- 
formed at the same time. 

As in the cases before, the Hadamard- and the tt/2- 
phase gate are performed only modulo a subsequent 
byproduct operator which is determined by the random 
measurement outcomes Sk 



Before we explain the functioning of the above 
gates, we would like to address the following questions: 
First, "How does one manage to occupy only those lattice 
sites with cluster qubits that are required for a particu- 
lar circuit but leaves the remaining ones empty?" . The 
answer to this question is that redundant qubits will not 
have to be removed physically. It is sufficient to measure 
each of them in the cr z -eigenbasis, as will be described in 
Section HTCl 

Second, "How can the described procedures for gate 
simulation be concatenated such that they represent a 
measurement based simulation of an entire circuit?" . It 
seems at first sight that the described building blocks 
would only lead to a computational scheme consisting 
of repeated steps of entangling operations and measure- 
ments. This is not the case. As will be shown in Sec- 
tion lll Dl the three procedures stated are precisely of such 
a form that the described measurement-based scheme of 
quantum computation can be decomposed into them. 

The third question is: "How does one deal with the 
randomness of the measurement results that leads to the 
byproduct operators i|2"4"fl. (jSUJ) and (|3"TJ) ?" . The appear- 
ance of byproduct operators may suggest that there is a 
need for local correction operations to counteract these 
unwanted extra operators. However, there is neither a 
possibility for such counter rotations within the described 
model of quantum computation, nor is there a need. The 
scheme works with unit efficiency despite the random- 
ness of the individual measurement results, as will be 
discussed in Section ITTeI 



C. Removing the redundant cluster qubits 

A cluster state on a two-dimensional cluster of rectan- 
gular shape, say, is a resource that allows for any com- 
putation that fits on the cluster. If one realizes a certain 
quantum circuit on this cluster state, there will always 
be qubits on the cluster which are not needed for its re- 
alization. Such cluster qubits we call redundant for this 
particular circuit. 

In the description of the QC C as a quantum logic net- 
work, the first step of each computation will be to re- 
move these redundant cluster qubits. Fortunately, the 
situation is not such that we have to remove the qubits 
(or, more precisely, the carriers of the qubits) physically 
from the lattice. To make them ineffective to the real- 
ized circuit, it suffices to measure each of them in the 
cr z -eigenbasis. In this way, one is left with an entan- 
gled quantum state on the cluster Cn of the unmeasured 
qubits and a product state on C\Cn, 

\^}) c ^\Z)c\c n ^\^W}) Cn , (32) 

with \Z) C \ Cn = (® ie c\c N \ s i)i,z) and s i tnc results 
of the <7 z -measuremcnts. The resulting entangled state 
\4>{k,'}) Cn on the sub-cluster Cn is again a cluster state 
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obeying the set of equations JIJ. This can be seen as 
follows. First, by definition we have 

® i+{ ~i ytai - W } ) g . 

(33) 

Using the eigenvalue equations l|T)l. we now insert a cor- 
relation operator with a G Cn into the r.h.s of i|33|) 
between the projector and the state, and obtain 

\Z) C \c N ® I0{ K '}) Cn = (-l)<A"< a >|Z) C \ Cjv ® |0 W} ) Cjv , 

(34) 

with the correlation operators 



K '(a) = (a) 



r(<0 



(35) 



c6nbgh(a)nC w 

and the set specifying the eigenvalues 



6enbgh(o)n(C\Cjv) 



s;, mod 2. 



(36) 



As the new correlation operators K'^ in l|34|) only act 
on the cluster qubits in Cn, the states \4>t K n) r again 
obey eigenvalue equations of type Q), i.e. 

^ /(o) l^{«'}) Cw = (-1)"'I*K}> C „. V a e Cat. (37) 

There are |Cjy| such eigenvalue equations for a state of 
ICjvI qubits. Thus, the state \(j>t K >\} r is specified by (pT7f> 
up to a global phase. 

From (|36J) we find that the redundant qubits have some 
remaining influence on the process of computation. After 
they have been measured, the random measurement re- 
sults enter into the eigenvalues that specify the residual 
cluster state |0{k'}) Cjv on the cluster Cn- However, any 
cluster state \4>{k'}) c is equally good for computation 
as stated in Section Til Al From i|14[l it follows that 



aeCr. 



(38) 



The Pauli phase flip operators that appear on the r.h.s. 
of equation i|3l3|) may be absorbed into the subsequent 
measurements. This allows us to adopt the following two 
rules in the further discussion 

1. The redundant cluster qubits are discarded. 

We only consider the sub-cluster Cn ■ f*$g\ 

2. We assume that n' a = for all a G Cn- 

This reduction will make a number of expressions such as 
those for the byproduct operators more transparent and 
it will also simplify the remaining part of the universality 
proof. 



D. Concatenation of gate simulations 

A quantum circuit on the QC C is a spatial and tempo- 
ral pattern of measurements on individual qubits which 
have previously been entangled to form a cluster state. 
To better understand its functioning we would like -as in 
the network model of quantum computation- to decom- 
pose the circuit into basic building blocks. These building 
blocks should be such that out of them any circuit can 
be assembled. In explaining the QC C in a network lan- 
guage, we can relate the building blocks of a quantum 
logic network -the quantum gates- to building blocks of 
QC c -circuits. 

The fact that quantum gates can be combined to quan- 
tum logic networks is obvious. But the statement that, 
for a QC c -computation, measurement patterns which 
simulate gates can simply be patched together to give 
the measurement pattern for the whole circuit requires a 
proof. This proof is given next. 

We begin by stating the general form of the procedures 
to realize gates and sub-circuits. The reason why these 
procedures work is explained in subsequent sections. To 
realize a gate g on the QC C consider a cluster C(g). This 
cluster has an input section Ci(g), a body Cm (g) and an 
output section Co{g), with 



C I (g)uC M (g)UC (g)=C(g) 
C I ( g )nC M (g) = 
C T (g)nC (g) = 
C M (g)nC (g) = 0. 



(40) 



The measurement bases of the qubits in Cm (<?), the body 
of the gate g 7 encode g. The general scheme for proce- 
dures to realize a gate g on a cluster C(g) is 

Scheme 1 Simulation of the gate g on C(g), acting on 
the input state |^)i n - 

1. Prepare the input state |'0in) on Ci(g) and the 
qubits in Cm(.9) U Co(ff) individually in the state 
1+) = |0)a; such that the quantum state of all qubits 
in C(g) becomes 

\*m)c(g) = \^in) Cl (g) ® (g) |+)fc- (41) 
keC M (g)uC (g) 

2. Entangle | v E'i 11 )c( 9 ) by the interaction 

S (C{g)) = gab^ ( 42 ) 

o,6SC(g)| b-ae-fd 

such that the resulting quantum state is \^ e )c N = 

^ (C(9)) |*in) C ( 9 ). 

3. Measure the cluster qubits in Ci(g) U Cm (5), 
i.e. choose measurement bases specified by £ 
S 2 , k G C/ (g) U Cm (flO and obtain the random mea- 
surement results Sfc such that the projector 



p(Cjr( S )UC M (s)) _ 



1 + (-l) s V fc • a 



(fc) 



(43) 



keCi(g)uc M (g) 
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is applied. The resulting state is \^ u.t)c N = 

p(Ci(g)UC u (g))\^\ c _ 

Putting all three steps of SchemeHtogether, the relation 
between |*i n )cjv and |^ , ut>c lv is 

|*out) C „ = P^^Cuig)) S (C(g)) |^. n)cN _ (44) 

As we will show later, the state | v I'out)c N has the form 



|*out)Cjv 



\sk)k,f k ® |^out}c (a)> ( 45 ) 

V feeCj(g)UC A /( S ) / 



where \sk)k,r k denotes the state of the qubit k after the 
observable f% ■ <f^ has been measured and the measure- 
ment outcome was s&, and 



l^OUt) = Uz,gUg\lpi n }. 



(46) 



Therein, U g is the desired unitary operation and C/s,g 
an extra multi-local rotation that depends on the mea- 
surement results {sfc | k £ Ci(g) U Cm (<?)}• The extra 
rotations Us , g are always in the Pauli group, i.e. 



(•*)-(■*)' 



(47) 



modulo a possible global phase, and n = \I\ = \0\. In 
(03 the cr 1 ' 1 denote Pauli operators acting on the logi- 
cal qubit i, not cluster qubit. The values Xi,Zi € {0, 1} 
are computed from the measurement outcomes {sfc | k £ 
C I (g)UC M (g)}. 

We now have all prerequisites at hand to explain why 
measurement patterns of basic gates can be combined to 
form the measurement pattern of the whole circuit which 
is obtained by combining the gates. The scheme of quan- 
tum computation with the QC C consists of a single entan- 
gling operation which creates the resource cluster state 
and, subsequently, of a scries of onc-qubit measurements 
on that state. We want to view the measurement pattern 
of a quantum circuit as being composed of basic blocks 
from whose function the function of the whole circuit can 
be deduced. To do so, we will explain computation on 
the QC C as a sequential process of performing the cir- 
cuit gate by gate. Then we have to demonstrate that 
the computational scheme as it is practically carried out, 
i.e. entangle once and afterwards only measure, and the 
sequential scheme that we use to explain the functioning 
of the circuit are mathematically equivalent. 

The sequential scheme is this. Consider a circuit 
U = ni^i U 9i that consists of a succession of gates 
5ii ■■i9W\ € N applied to some input state IV'in), leading 
to an output state |V>out) = U\^p- m ) which is then mea- 
sured. N denotes the network, i.e. the set of gates plus 
a description of their relation. For simplicity let us first 
assume that each gate gi acts on all of the logical qubits. 
Subsequently we will drop this assumption. 
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FIG. 3: Here the exchange of the order of the measurements 
and the entanglement operations is shown. Fhe crosses "x" 
denote the one-qubit measurements and the horizontal lines 
between adjacent cluster qubits denote the unitary transfor- 
mations S a ' a+1 . 



First, a quantum state \ip- m )i ® |+)cjAi ^ s prepared. 
Then each gate, one after another, is realized on a sub- 
cluster C(g) C Cn according to Scheme H Finally, the 
output is measured as usual. 

In the step of carrying out the gate gi the state of the 
quantum register is, besides being processed, also tele- 
ported from Ci(gi) = C (gi-i) to C {gi) = Cj(g i+ i). In 
this way, by carrying out the gate gi the input for the 
successor gate gi+\ is provided. To proceed with the re- 
alization of gate gi+i, in accordance with Scheme ^ the 
sub-cluster C(gi+\) is entangled via SCCfe+i^and subse- 
quently the cluster qubits in C/(_g,+i)UCM(s'i+i) are mea- 
sured. This completes the realization of gate <jj+i and at 
the same time writes the input for gi+2, and so on. 

The reason why the sequential scheme just described 
is equivalent to the entangle-once-and-then-only-measure 
scheme is the following. The entanglement operations at 
the various stages of the sequential scheme commute with 
all the measurements carried out earlier. This holds be- 
cause both operations, entangling operation and earlier 
measurement, act on different particles. Thus, the oper- 
ations may be reordered in such a way that in a first step 
all entangling operations S^- C ^ 9i ^ act on the initial state 
and afterwards all the measurements are performed. 

The exchange of the order of the one-particle measure- 
ments and the two-particle Ising interactions is shown 
in Fig. |21 for a ID cluster. In one dimension the de- 
composition of a cluster into sub-clusters, as displayed 
in Fig. |3| is clear. However, the interesting cases for 
QC c -computations are clusters in 2D and 3D; and there 
we must state more precisely what "decomposition of a 
cluster into sub-clusters" means. The use of basic notions 
from graph theory will prove helpful for this purpose. 

To decompose a cluster into sub-clusters means in 
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more precise terms to decompose the associated graph 
G(Cn, Ec N ) into subgraphs. That is, we have to decom- 
pose both the vertices and the edges of the graph. Each 
vertex a G Cn has to belong to a subset C(gi), where 



--JV 



(48) 



and the sets C{gi) of vertices corresponding to the gates 
<?i may overlap on their input- and output vertices. 

Correspondingly, the set Eq n of edges, defined in the 
same way as Eq in J22J), is decomposed into subsets 



E c 



\M\ 



(49) 



but the subsets E{g^ of edges are not allowed to overlap, 



Vi,i = 1. J = E( 9i )nE( 9j ) = 



(50) 



The rules for the decomposition of edges (|49fl and i|50|) 
are, as we shall see, central for the universality proof. 

Further, for the decomposition to be useful, the subsets 
C(gi) and P(<?i) must fulfill a number of constraints. The 
first of these is that each pair (C(gi), E(gi)) is again a 
graph, G(C(gi), E(tji)). This requires, in particular, that 
the endpoints of all the edges in E{gi) are in C{gi), 



Vi = l,..,\M\:E(g i )cC(g i )xC(g i ). 



(51) 



For details on the graph decomposition, in particular for 
conditions on the subgraphs imposed to guarantee l|49|) 
and H5Ufl sec Appendix [XJ 

Now consider the concatenation g 2 °gi of the two gates 
<7i, realized on a cluster C{g\), and gi-, realized on a 
cluster C(<72), each of them by a procedure according to 
scheme ^ The composite circuit g = 92 ° gi is realized 
on the cluster 



with 



C(g)=C( gi )UC(g 2 ), 



Cj(g) = Cjfoi) U (C/(<? 2 )\Co(<7i)) 
Co(g) = CoGte) u (CoteOVMsa)). 



(52) 



(53) 



Now, the procedure to perform the two gates <?i, 32 
sequentially is 

1. Prepare the state |*i n )c( s ) = IVWcjfo) ® 
l+)c( 9 )\c,(s)- 

2. Entangle the qubits on the sub-cluster C{gi) via 

Si := 5 (c(si)) = Yl S ab . (54) 

a,b£C(gi)\ b-ae^td 



3. Measure the qubits in Ci(gi) UCa/(<7i), resulting in 
the projector P Cl ( gi )uc M (g 1 ) =■ Pi, 



P = 



1 + (-l) st r fe • 



(55) 



fceC/( ff i)UC M (gi) 



Therein is the outcome of the measurement of 
qubit k and r% the respective measurement direc- 
tion. 



4. Entangle the qubits on the sub-cluster C(g2) via 

2b (56) 



S 2 := S< c ^» = 5 a 

a,bGC(g 2 )| b-aEld 



5. Measure the qubits in £7(52) U Cm (52), resulting in 
the projector P C ;(s2)uC„( ff2 ) =: P 2 , 



Po 



1 + (-l) Sfc 7^ - a 



fc6C 7 (g 2 )UC Jf (g 2 ) 

The procedure results in an output state 



I* 



out )C(g) 



P2S 2 P 1 S 1 \t> i 



(57) 



(58) 



that has the form 



I* 



out/C(ff) 



\sk)k,r k ® l^out)co(s)' ( 59 ) 



,C/( S )uCm( 9 ) 



with 



(60) 



according to l|46|) . 

As will be shown next, the above procedure is equiv- 
alent to a procedure of Scheme ^ applied to the cluster 
C(g) = C(gi) U C{g2), i.e. when, first, all qubits in C(g) 
are entangled and, second, all but the output qubits of 
Co(g) are measured. 

The procedure according to Scheme ^ yields the state 

\Kut)c {g) = P^ uc ^» S^» |* in ) c(g)l (61) 

and we now have to show that the output states 
|*out)c(g) m <EH1) and | v I / out)c(5) m ijSTJl are the same for 
all input states \^i n }c( g )- 

First note that the operations Pi and S2 commute 
since they act on different particles. P\ acts on the qubits 
in Cf(<?i) U Cm(<?i) while ^2 acts on C(<?2)- The sub- 
clusters associated with the gates may overlap only via 
their input- and output qubits. This is intuitively clear, 
and also follows from the decomposition constraint 1|A5() . 
As the gate g\ is applied before 172, of C(<?i) only the 
qubits in Co{gi) may overlap with the qubits in Cj^)- 
Thus, (Ci(g 1 )UCu{gi))nC(g 2 ) = 0. Therefore 



P2 S2 Pi S\ — P2 Pi S2 Si . 



(62) 
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Now note that as a direct consequence of ((53(1 the union 
of the input- and body section of the composite gate g on 
the cluster C(g) are made up by the union of the input- 
and body sections of the two individual gates g\ and g 2 , 
i.e. 

Ci(gi)L>C M {g 1 )UC I (g 2 )UC M (g 2 )=C I (g)uC M (g). (63) 

Further, from the decomposition constraint l(A5(l and 
from the fact that g\ is applied before g 2 it follows that 
the input- and body sections of gates g\ and g 2 do not 
intersect, 

(C/Ga) U C M {gi)) n (d(g 2 ) U C M {g 2 )) = 0. (64) 



Therefore, 
P 2 P 



1 + (-l) Sk f k ■ ffW 



Cj(s 2 )UC M (S2) 



(fc) 



(65) 



Cj( S1 )UCm(si) 

^ 1 + (-l) a -r k ■ 

Cj( ff )UCAj(£j) 

_ p(Ci{g)uc M (g)) 

where the second line holds by 1(63(1 and (|64|l . We find 
that measurement patterns corresponding to the projec- 
tions Pi and P2 can be patched together to form the 
measurement pattern on the cluster C(g). 

The same holds for the entangling operations. The en- 
tangling operation Si on C(<?i) and S 2 on C(g 2 ) combined 
give the entangling operation S^^ 9 " on C(g), 

S 2 S 1 = S^ g) \ (66) 

because of the central rule (|5()ll . 

Inserting ((65(1 and 1)66(1 into ((62(1 yields 

P 2 S 2 Pi S*i = p( c '(3)uc M (g)) s (C(g))^ ( 67 ) 

and therefore, if we compare 1(58(1 and 1(61(1 we find that 
|5 , out )c( ff ), the output state of the sequential realization 
of the two gates g\ and g 2 , and | v I / out)c(g) i the output 
state of the standard procedure applied to the composite 
circuit, are indeed the same for all inputs |^ / i n )c(g)- Thus 
both realizations, the sequential and the non-sequential, 
are equivalent. 

This composition can be iterated so that the entire 
circuit can be realized via the standard procedure of 
Scheme ^ The measurement pattern of the circuit is 
thereby obtained by patching together the measurement 
patterns of the gates the circuit is composed of. 

From ((60(1 it follows that the quantum input \ipi n ) and 
the quantum output \ipout) of the unitary evolution are 
related via 



\lpo 



(68) 



The random but known byproduct operators Us t9i that 
appear in 1(68(1 are dealt with in Section 111 El The gates 
gi € A/" are labeled corresponding to the order of their 
action. 

Now, we want to specify to the case where the quan- 
tum input is known and where the quantum output is 
measured. This is the situation which interests us most 
in this paper. Examples of such a situation are Shor's 
factoring algorithm and Grover's search algorithm. In 
both cases, the quantum input is \tpi n ) = ®j=i 

Let us denote the input section of the whole cluster 
C, comprising the input qubits of the network simula- 
tion, as /; and the output section, comprising the qubits 
of the readout quantum register, as O. As long as the 
quantum input is known it is sufficient to consider the 
state |+)/ = (& ieI F° r different but known input 

states \ipin)i one can always find a transformation U{ n 
such that \ip- m )i = U- m \+)i and instead of realizing some 
unitary transformation U on \ijjm)i one realizes UU- m on 
1+)/- ' 

Preparing an input state |+)/ and entangling it 
via is the same as creating a cluster state \4>)c, 

|+)/ ® \+) C \i = S (c) |+)c = \4>)c ■ This holds be- 
cause the state S^ c ' \+}c obeys the eigenvalue equations 
(Q and, as we have stated earlier, these eigenvalue equa- 
tions determine the state completely. Thus the created 
state is a cluster state \<p)c that could as well have been 
prepared by any other means. 

Once the quantum output is read then all clus- 
ter qubits have been measured. Therefore, the entire 
procedure of realizing a quantum computation on the 
QC C amounts to 

Scheme 2 Performing a computation on the QC C . 

1. Prepare a cluster state \<fi{ K y) c of sufficient size. 

2. Perform a sequence of measurements on |0{ K }) C 
and obtain the result of the computation from all 
the measurement outcomes. 

The link between the network model and the QC C is 
established by Scheme QJ The elementary constituents of 
the quantum logic network are mapped onto the corre- 
sponding basic blocks of the QC C . In this way, Scheme^ 
helps to understand why the QC C works. 

However, a QC c -computation is more appropriately 
described by Scheme |21 than by Scheme Scheme [5] 
does not use the notion of quantum gates, but only of 
a spatial and temporal measurement pattern. Once uni- 
versality of the QC C is established, to demonstrate the 
functioning of specific QC C -algorithms one would prefer 
decomposing a measurement pattern directly into sub- 
patterns rather than decomposing a network simulation 
into simulations of gates. A tool for the direct approach 
is provided by Theorem 2] in Section III Fl 
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E. Randomness of the measurement results 

We will now show that the described scheme of quan- 
tum computation with the QC C works with unit effi- 
ciency despite the randomness of the individual measure- 
ment results. 

First note that a byproduct operator Us that acts af- 
ter the final unitary gate U 9]N] does not jeopardize the 
scheme. Its only effect is that the results of the readout 
measurements have to be reinterpreted. The byproduct 
operator Us that acts upon the logical output qubits 1 . . n 
has the form 



Us = U 



(69) 



where x%,Zi € {0, 1} for 1 < i < n. Let the qubits on the 
cluster which arc left unmeasured be labeled in the same 
way as the readout qubits of the quantum logic network. 

The qubits on the cluster which take the role of the 
readout qubits are, at this point, in a state t/s|out), 
where |out) is the output state of the correspond- 
ing quantum logic network. The computation is com- 
pleted by measuring each qubit in the cr z -eigenbasis, 
thereby obtaining the measurement results {s^}, say. 
In the QC c -schcmc, one measures the state J7s|out) 
directly whereby outcomes {s^} arc obtained and the 
readout qubits arc projected into the state \M) = 



n 



n l+(-l) s 



E/s|out). Depending on the byproduct 



1 2 

operator Us, the set of measurement results {s} in gen- 
eral has a different interpretation from what the network 
readout {s^} would have. The measurement basis is the 
same. From (|fi9|l one obtains 



\M) 



n ( " } - w 



A rr ! + (-!) °Y 




Us |out) (70) 



out) 



(=1 



From (|70|) we see that a <7 z -measurement on the state 
Us |out) with results {s} represents the same algorithmic 
output as a cr 2 -mcasurement of the state |out) with the 
results {s'j}, where the sets {s} and {s£} are related by 



3d 2. 



(71) 



The set {s^} represents the result of the computation. 
It can be calculated from the results {si} of the a z - 
mcasurcments on the "readout" cluster qubits, and the 
values {xi} which are determined by the byproduct op- 
erator Us- 

Thus we find that one can cope with the randomness of 
the measurement results provided the byproduct opera- 
tors Us, gi in CIH can be propagated forward through the 



subsequent gates such that they act on the cluster qubits 
representing the output register. This can be done. To 
propagate the byproduct operators we use the propaga- 
tion relations 



CNOT(c, t)atP 
CNOT(c, t)<j { x c) 
CNOT(c, t)ai t] 
CNOT(c, t)a z c) 

for the CNOT gate, 

Unot[£, V,(}<7x 



f ri t) CNOT(c, t) 
ai c) ai t) CNOT(c, t) 
CTi c) f ji t) CNOT(c,i) 

CTi c) CNOT(c, t) 



<rxU Ro t[(,,-ri,C} 



(72) 



(73) 



VzU Rot [-£,ri,-(] 
for general rotations U Ro t[£,, rj, £] as defined in (|25|) . and 



Ha x 
Ha z 
U z [ir/2}a x 
U z [ir/2]a z 



a z H 

<j y U z [ir/2] 
cr z U z [n/2] 



(74) 



for the Hadamard- and 7r/2-phase gate. The propaga- 
tion relations i|73|) apply to general rotations realized via 
Procedures-including Hadamard- and 7r/2-phase gates- 
while the propagation relations (|74|> apply to Hadamard- 
and 7r/2-phase gates as realized via Procedure 

Note that the propagation relations l|72(l - (|74(l are such 
that Pauli operators are mapped onto Pauli operators 
under propagation and thus the byproduct operators re- 
main in the Pauli group when being propagated. Fur- 
ther note that there is a difference between the relations 
for propagation through gates which are in the Clifford 
group and through those which are not. For CNOT-, 
Hadamard- and 7r/2-phase gates the byproduct operator 
changes under propagation while the gate remains un- 
changed. This holds for all gates in the Clifford group, 
because the propagation relations for Clifford gates are 
of the form U g Us = {UgUsU^Ug as J2} and dJ», i.e. 
the byproduct operator Us is conjugated under the gate, 
and the Clifford group by its definition as the normalizcr 
of the Pauli group maps Pauli operators onto Pauli oper- 
ators under conjugation. For gates which are not in the 
Clifford group this would in general not work and there- 
fore, for rotations which are not in the Clifford group, 
the propagation relations are different. There, the gate 
is conjugated under the byproduct operator; and thus 
the byproduct operator remains unchanged in propaga- 
tion while the gate is modified. In both cases, the forward 
propagation leaves the byproduct operators in the Pauli 
group. In particular, their tensor product structure is 
maintained. 

Let us now discuss how byproduct operator propaga- 
tion affects the scheme of computation with the QC C . In 
Section III Dl we arrived at the conclusion 168|) that by 
patching the measurement patterns of individual gates 
together and keeping the measurement bases fixed, we 
can realize a composite unitary evolution on some input 



state \tpin), 



out/ 



where f/ ffi is the i-th unitary gate in the circuit and U^ gi 
the byproduct operator resulting from the realization of 
that gate. Now using the above propagation relations, 
(|68[1 can be rewritten in the following way 



out) 




(75) 



Therein, U-£, gi \ n are forward propagated byproduct op- 
erators resulting from the byproduct operators Us. gi of 
the gates gi. They accumulate to the total byproduct op- 
erator Uy, whose effect on the result of the computation 
is contained in (|71|) . 



IN 
O- 



u. 



u, 



12 

OUT 



"1 



U, 







1 


b 







U x --Q 



U 



o, 



FIG. 4: Vertical cuts. The vertical cuts intersect each qubit 
line exactly once but do not intersect gates. Thus, Oi, Oj and 
Q are vertical cuts, but is not. The cut Oi intersects the 
rotation U x just before its input. For two of the rotations in 
the displayed network, the sub-clusters on which these gates 
are realized are symbolically displayed in gray underlay. Via 
the measurement of the cluster qubits a and 6 (displayed as 
black dots with white border), the rotation angles of the re- 
spective rotations U x and U z are set. 



C^ = n^.«ln- ( 76 ) 

i=i 

Further, the U' g . are the gates modified under the prop- 
agation of the byproduct operators. As discussed above, 
for gates in the Clifford group we have 

U' g = U g , V g e Clifford group, (77) 

as can be seen from (17211 and (|74|) . 

Gates which are not in the Clifford group are modi- 
fied by byproduct operator propagation. Specifically, the 
general rotations i|25|) arc conjugated as can be seen from 
(I73|l . From the structure of (|68fl we see that only the 
byproduct operators of gates gk earlier than gi in the net- 
work may have an effect on U gi , i.e. those with k < i. To 
give an explicit expression, let us define | q , which 

are byproduct operators J7s ;Sfc propagated forward by the 
propagation relations (|72"|) - l(74"|) to the vertical cut Oi 
through the network, see Fig.0] A vertical cut through a 
network is a cut which intersects each qubit line exactly 
once and does not intersect gates. The vertical cut Oi 
has the additional property that it intersects the network 
just before the input of gate gi. The relation between a 
rotation U' g . modified by the byproduct operators and 
the non-modified rotation U gi is 

k = fn u ^\o\uJ n u s , gk \ \ , 

\fe|fe<-i / \k\k<i J (78) 

V U Bi G SU(2). 

Now that we have investigated the effect of byproduct 
operator propagation on the individual gates let us re- 
turn to equation (|75f) . There, we find that the opera- 
tions which act on the input state \ipin) group into two 



factors. The first is composed of the modified gate op- 
erations U'g. and the second of the forward propagated 
byproduct operators. The second factor gives the accu- 
mulated byproduct operator U% and is absorbed into the 
result of the computation via l|71|) . It does not cause any 
complication. 

So what remains is the first factor, and we find that the 
unitary evolution of the input state iV'in) that is realized 
is composed of the modified gates U' g . . The gates we will 
realize are thus the U', not [/„.. However, the standard 
procedures 1 - 3 in Section III Bl are for the operations 
U gi . Thus we have to read l|78|l in reverse. We need to 
deduce U 9i from U' . Once the gates gk for all k < i have 
been realized, this can be done for each gate gi since the 
byproduct operators are then known for all k < i. 
Finally, with U 9i determined from U' g . , Procedure Ogives 
the measurement bases required for the realization of the 
gate gi. Please note that it is a sufficient criterion for 
the realization of the gate gi that all gates gk with k < i 
must have been realized before, but not a necessary one. 

Let us, at this point, address the question of tempo- 
ral ordering more explicitly. For proper discussion of 
the temporal ordering we have to step out of the net- 
work frame for a moment. First note that in case of the 
QC C the basic primitive are measurements. Thus, the 
temporal complexity will be determined by the temporal 
ordering of these measurements, unlike in quantum logic 
networks, where it depends on the ordering of gates. The 
most efficient ordering of measurements that simulates a 
quantum logic network is not pre-described by the tem- 
poral ordering of the gates in this network. 

A temporal ordering among the measurements is in- 
ferred from the requirement to keep the computation on 
the QC C deterministic in spite of of the randomness in- 
troduced by the measurements. This randomness is ac- 
counted for by the byproduct operators. The key to ob- 
tain the temporal ordering of measurements is eq. I|78|l . 
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There, the byproduct operators U^. gk \ o i may modify Eu- 
ler angles of the one-qubit rotations in the network and 
consequently change measurement bases. The tempo- 
ral ordering thus arises due to the fact that bases for 
one-qubit measurements must be chosen in accordance 
with outcomes obtained from the measurements of other 
qubits. 

For each cluster qubit q that needs to be measured 
in a non-trivial basis, i.e. not in the eigenbasis of <j Xl 
o y or a zi a set of cluster qubits pi can be identified, 
whose measurement outcomes influence the choice of the 
measurement basis for qubit q. We say that q is in the 
forward cone Q of pi, q £ k(pi). Each cluster qubit has 
a forward cone, and in no forward cone there appears a 
qubit which is measured in a trivial basis. 

The rule is that a cluster qubit q can only be mea- 
sured once all cluster qubits pi for which q £ fc(pi) have 
been measured earlier. The forward cones thereby gener- 
ate an anti-reflexive partial ordering among the measure- 
ments from which the most efficient measurement strat- 
egy can be inferred, see Q. Gates in the Clifford group 
do not contribute to the temporal complexity of a QC C - 
algorithm, see Section UTTl 

F. Using quantum correlations for quantum 
computation 

In this section we give a criterion which allows do 
demonstrate the functioning of the QC c -simulations of 
unitary gates in a compact way. 

Before we state the theorem, let us make the notion 
of a measurement pattern more precise. In a QC C - 
computation one can only choose the measurement bases, 
while the measurement outcomes are random. This is 
sufficient for deterministic computation. Thus one can 
perform measurements specified by a spatial and tempo- 
ral pattern of measurement bases but one cannot control 
into which of the two eigenstates the qubits are projected. 

Definition 1 A measurement pattern on a cluster 

C is a set of vectors 

M {c) = {f a eS 2 \aeC}, (79) 

defining the measurement bases of the one-qubit measure- 
ments on C. 

If this pattern M.^ of measurements is applied on an 
initial state \*&£ )c and thereby the set of measurement 
outcomes 

{s}c = {s a S {0, 1} | a £ C} (80) 

is obtained, then the resulting state | v I'ai)c is, modulo 
norm factor, given by \~$m)c = PijUM) l*£>c, where 

'ffCM)-(8 1 + ( - 1) :*-' m . («) 



Besides, let us introduce some conventions for labeling. 
BeCj(g) and Co{g) such that |Cr(g)| = |Co(.g)| = « where 
n is the number of logical qubits processed by g. Opera- 
tors acting on qubits p <G Ci(g) and q £ Co{g) are labeled 
by upper indices {Ci{g),i) and {Co{g),i'), 1 < i, i' < n, 
respectively. The qubits p € Ci(g) and q £ Co(g) are 
ordered from 1 to n in the same way as the logical qubits 
that they represent. 

We make a distinction between the gate g and the 
unitary transformation U it realizes. The gate g £ Af 
does, besides specifying the unitary transformation U, 
also comprise the information about the location of the 
gate within the network. 

After these definitions and conventions we can now 
state the following theorem 

Theorem 1 Be C(g) = d{g) U C M (g) U C (g) with 

d( g ) n C M {g) = dig) n C (g) = C M (g) n C (g) = 

a cluster for the simulation of a gate g, realizing the uni- 
tary transformation U , and \4>)c( g ) the cluster state on 
the cluster C(g). 

Suppose, the state |V)c( s ) = P^'^iM) \4>)c{ a ) obeys 
the 2n eigenvalue equations 

(^%t) (Co(s) V) C(ff) = (-1) A -I^)c( 9 ), 

(82) 

with \ Xy i, X z .i £ {0, 1} and 1 < i < n. 

Then, on the cluster C(g) the gate g acting on an arbi- 
trary quantum input state \ipi n } can be realized according 
to Scheme^\with the measurement directions in Cm(s) 
described by M^ Cm ^ 9 ^ and the measurements of the qubits 
inCi(g) being a x -measurements. Thereby, the input- and 
output state in the simulation of g are related via 

l^out) = UUx hM, (83) 
where is a byproduct operator given by 

n 

U^= (g) (al^+^aW)^ 1 - (84) 

(Ci{g)3i)=l 

The significance of the above theorem is that it pro- 
vides a comparatively simple criterion for the functioning 
of gate simulations on the QC C . 

In Scheme^ after read- in of the input state and the en- 
tangling operation S 1 -- ^, i.e. before the measurements 
that realize the gate are performed, the resulting state 
carries the quantum input in an encoded form. This state 
is in general not a cluster state. It is therefore not clear a 
priori that cluster state correlations alone are sufficient to 
explain the functioning of the gate. However, this is what 
Theorem n states. To prove the functioning of a gate g 
realized via Scheme ^ it is sufficient to demonstrate that 
a cluster state on C(g) exhibits certain quantum corre- 
lations. About the variable input one does not need to 
worry. 
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This is convenient in two ways. First, we can base the 
explanation of the gates directly on the eigenvalue equa- 
tions which were also used to define the cluster states 
in a compact way. The quantum correlations required to 
explain the functioning of the gates are derived from the 
basic correlations J5J) rather easily and thus the use of 
Theorem n makes the explanation of the gates compact. 

Second, Theorem ^ is a tool to demonstrate the func- 
tioning of QC c -circuits without having to repeat the 
whole universality proof for each particular circuit un- 
der consideration. Scheme [5] describes the computation 
as a series of one-qubit measurements on a cluster state. 
An accordance with this, instead of decomposing a cir- 
cuit simulation into gate simulations as done in Scheme^ 
a measurement pattern is decomposed into sub-patterns. 
The effect of these measurement sub-patterns is tested 
via the criterion (|82|l in Theorem ^ 

Before we turn to the proof of Theorem let us note 
that the measurements described by P^\ M (-M(g)), 
as they have full rank, project the initial cluster 
state \<p)c( g ) into a tensor product state, \ip)c(g) = 
l m )cjv,(s) ® w)ci(g)uCo(g)- Thereof only the second fac- 



tor, 



/Cj(g)UCo(s)' 



is of interest. This state alone sat- 



isfies the eigenvalue equations 1|82|) . and is uniquely de- 
termined by these equations. To see this, consider the 

state W)ci(g)uc (g) = t/t K'>c / (s)uc (9)- :t satisfies the 
2n eigenvalue equations 



(85) 



where we have written in short for \^')ci(g)uCo{g)- 
The state \^')ci(g)uc {g) i s uniquely defined by the above 
set of commuting obscrvables, it is a product of Bell 
states. Therefore, \ip)ci( g )uc (g) IS uniquely defined as 
well. 

Proof of Theorem QJ We will discuss the functioning 
of the gates for two cases of inputs. First, for all input 
states in the computational basis. This leaves relative 
phases open which have to be determined. To fix them, 
we discuss second the input state with all qubits indi- 
vidually in |+). As we will sec, from these two cases it 
can be concluded that the gate simulation works for all 
input states of the computational basis. This is suffi- 
cient because of the linearity of the applied operations; if 
the gate simulations work for states of the computational 
basis then they work for superpositions of such inputs as 
well. 

Case 1: The input \ip- m ) is one of the states of the 
computational basis, i.e. \ipi n ) = |z) := <SC=i \%i)z,i with 
Zi G {0, 1}, i = 1 .. n. Then the state \^ out (z))c( g ) of the 
qubits in C [after performing a procedure according to 
Scheme^ using a measurement pattern M.(° M (s)) on the 
body Cm(9) of the gate g, and applying cr^-measurements 



on Ci(g)} is 
n (z) |* out (z)) c(ff) 



P 



f s ^\ X )pl c s ;'^M)s^ \z) Cl(g ^\+) CMig)UCoia) , 

(86) 

with norm factors no(z) that are nonzero for all z, as we 
shall show later. 

The input |z) in l|86|) satisfies the equation 



m(z)\z)=P£ 



(C(g)) 



(87) 



with P, 



(Ci(ff)) 



<8>IU 



, and nj(z) = 1/2 



n/2 



for all z. Now note that 5 (c(s)) and Pjf z l(s)) , as well 



,(Cm(9)) ( 



, {s} >(M) and Pf^ ia) \ commute. Thus, |* out (z)) c(g) 
can be written as 



n' (z) |*out(z))c(s) = 



P, 



(Ci(ff)) 



W 



P, 



(C/(ff)) 



{X)P : 

POP 



{dig)) p (C„(s)) 



(Cj(fl)) 



C(fl)> 



where 



is specified by the eigenvalue equations jBSj) 



in Theorem ^ 

Let us, at this point, emphasize that the projections 
P { ( f } /(s)) (A) and Pf^ g)) in (gg) are of very different 

origin. The projector Pj?y 9 ''(X) describes the action 
of the a x -measurements on the qubits in Ci(g). These 
measurements are part of the procedure to realize some 
gate g on the cluster C(g). One has no control over the 
thereby obtained measurement outcomes {s} specifying 



P { ¥ s)) (l). In contrast, the projector P^ 9 " docs not 
correspond to measurements that are performed in real- 
ity. Instead, it is introduced as an auxiliary construction 
that allows one to relate the processing of quantum in- 
puts to quantum correlations in cluster states. The pa- 
rameters z specifying the quantum input |z) and thus the 

projector Pf^ in JH3 can be chosen freely. 

The goal is to find for the state ont {z)) c( g ) an ex- 
pression involving the transformation U acting on the 
input |z). To accomplish this, first observe that for the 
state on the r.h.s of l|88|) via H82|) the following eigenvalue 
equations hold 



(Ci(g)) 



(Co) 



,(Ci(g)), 



W 



pf4 i9) \x)p : 



(C,(g)) 
Z.z 



>C(g) 



(89) 
with i = l..n. 

To make use of the equations (|89|l we need to prove 
that P ( { C s} l{9) \x) P|^ (s)) \tp)c(g) ^ for all z under the 
assumptions of theorem ^ 



For this, we consider the scalar 



C(g)\ 



IP 



(C/ (<?)), 



Z.z 



)C(g) 
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and write P^z m th° f° rm 



P, 



(Ci(g)) _ 



Z.z 



(Ci(g)) 



(90) 



where C Ci ^ Vfc = 1..2". For each Ik we choose 
an i G and insert the respective eigenvalue equation 
from the upper line of into c(g){^\ <8> je / fc ^W)c{g)- 
Since ®,^<r« and (t/^Vt) (C ° (ff)) anti _ 

commute, c( ff )(^| <8> je / fc ^ IV>}e( ff ) = for all I k . Thus, 
with JU, one finds c(ff) (ip\P^ (g)) \ip)c(g) = 1/2", such 
that P^z |t/')c(s) 7^ and therefore also 

p(c,( 3 )) Wp ^( 9 ))| V;)c(ff) ^ 05 (91) 

or, in other words, n' c (z) ^ for all z. 

Due to the fact that the projections -P^z ana - 
■P/^ are of full rank the above state has the form 



p(Ci(g))(X\ p( c Hs)) 



n' (z) |s)xA(s) ® \ m )c M (g) ® IV'out(z)) Co (9)1 

(92) 

n 

where (s^c, = <S> \$i)x,h and l m )c Af ( 9 ) is some prod- 

(Ci3i)=l 

uct state with || |m)c M ( ff )|| = 1. Elaborating the argu- 
ment that leads to 1(91(1 one finds that n' (z) = 1/2™ and 
no(z) = 1/2' 1 / 2 , but at this point the precise values of 
the normalization factors are not important as long as 
they are nonzero. 

In (|92|) only the third factor of the state on the r.h.s. 
is interesting, and this factor is determined by the eigen- 
value equations l(89|): 



|Vw(z)> =e^ z >£/[/ E |z), 



(93) 



where Us is given by Now, because of (f%£|) with 

n' G (z) 7^ Vz, a solution (jHU with i|9"3|> for the state 
p(C/(»)) (X) P^ C z j(9)) |i/')c( g ) is also a solution for the state 
l^'out(z))c( s ) , and one finally obtains 

|*out(z))c(g) = e"' (z) |s) :!; ,c J (9)®l m )cM( s )®[C / C / s|z)] Co(g) . 

(94) 

There appear no additional norm factors in (|94|l because 
the states on the l.h.s. and the r.h.s. are both normalized 
to unity. 

The solution 1(94(1 still allows for one free parameter, 
the phase factor e"'( z \ Note that, a priori, the phase 
factors for different z can all be different. 

This concludes the discussion of case 1. We have found 
in l|94|) that the realized gate acts as 



where the gate D is diagonal in the computational basis 
and contains all the phases e ir, ( z >. What remains is to 
show that D = 1 modulo a possible global phase. 

Case 2. Now the same procedure is applied for the 
input state \ipi n ) = |+) := Then, the state 

| v I / out(+))c(g) that results from the gate simulation is 



no(+) |*a»t(+)) C(fl ) = Pl C 4 i9)) (X)P^(M) \4>) c{g) , 

(96) 

with a nonzero norm factor no(+). Using the upper line 
of eigenvalue equations JH2J), the state \^out(+))c(g) is 
found to obey the eigenvalue equations 

(Uo^uf° {9)) |*out( + )) C(s ) = {-lf- i+Si |*out(+)) C(fl ). 

(97) 

The eigenvalue equations (|97|l in combination with 196|) 
imply that 

|*out(+))c( 9 ) = e^ls),,^^ ® \m)c M{g) <8 [UUs\+)] Co(g) , 

(98) 

with x being a free parameter. Therefore, on the input 
state |+) the gate simulation acts as 



U = e lx U Us 



(99) 



This observation concludes the discussion of case 2. 

The fact that l|94|l and (|98|) hold simultaneously im- 
poses stringent conditions on the phases t](z). To see 
this, let us evaluate the scalar product 

c x = c(g)(^out{+)\UUs\s) x _ Clig) <g> \m) CM{g) <g> \+) Co (g)- 

(100) 

From (|98|l it follows immediately that 

(101) 

On the other hand, since |+) = 1/2' 1 / 2 J2 Z £{q u» l z ) and, 

by linearity, |* TOt (+)> = l/2" /2 E zG{ o,i}" l*out(*)>> 
from it follows that 



1 E 

2™ ' 

ze{0,l}" 



(102) 



The sum in i(102|> runs over 2" terms. Thus, with 
| e — t»j(a)| _ 2 f or a n Zj it follows from the triangle in- 
equality that \c x \ < 1. The modulus of c x can be unity 
only if all e _M) ( z ) are equal. As (|101fl shows, \c x \ is indeed 
equal to unity. Therefore, the phase factors e 171 ^ must 
all be the same, and with 11101(1 and 1(102(1 , 



(103) 



U = UU^D 



(95) 



If we now insert (|103|) into H94|) we find that the gate sim- 
ulation acts upon every input state in the computational 
basis, and thus upon every input state, as U g = e lx U Us- 
Therein, the global phase factor e %x has no effect. Thus 
we find that the gate simulation indeed acts as stated in 
JH3 and IS). □ 
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We would like to acknowledge that a similar theorem 
restricted to gates in the Clifford group has been obtained 
in P3. 

Let us conclude this section with some comments on 
how to use this theorem. First, note that Theorem^ does 
not imply anything about the temporal order of measure- 
ments within a gate simulation. In particular it should 
be understood that a procedure according to Scheme ^ is 
not such that first the measurements on the cluster qubits 
in Cm {9) an d thereafter the measurements in Ci{g) are 
performed. 

Instead, first all those cluster qubits q e Cj(g) UCm(<?) 
are measured whose measurement basis is the eigenbasis 
of either a x or a y (remember that, after the removal of 
the redundant cluster qubits as described in Section fll CI 
we are dealing with clusters Cm such that, apart from 
the readout, no measurements in the a z -eigenbasis oc- 
cur). Second, possibly in several subsequent rounds, the 
remaining measurements are performed in bases which 
are chosen according to previous measurement results. 

Let us now discuss how to choose the appropriate mea- 
surement bases. First note that the unitary operations 
Us and U in (|83|) both depend on measurement results 
of qubits in C(g), 

Us = Us({s i \i£C I (g)UC M (g)}), 
U = u(M^\{ Si \i€C M (9)}). (1 ° 4) 

The dependence of Us on {si\ i £ C},[{g)} arises through 
the {\ x .i, X z .i\i = 1 .. n} of if%5)) . 

Now note that in (j83|) the order of the unitary gate U 
and the byproduct operator Us is the opposite of what is 
required in (|68|l . Therefore, the order of these operators 
has to be interchanged, which is achieved by propagat- 
ing the byproduct operator Us through the gate U. For 
gates or sub-circuits given as a quantum logic network 
composed of CNOT-gates and one-qubit rotations, this 
task can be performed using the propagation relations 
C2J, O and HU). The result is 

U = U {M^^\{s} Cm }) Us ({s} Cl uc M ) = 

Us({s} Cl uc M )U' (M^\{s} CiUCm ). 

(105) 

Now, the choice of measurement bases in jV/( CM (f)) is 
allowed to be adaptive, that is the measurement bases 
may depend on measurement outcomes at other cluster 
qubits, M^ Cm{9)) = M (Cm{9)) {{s}). For the realization 
of the gate U g , the measurement bases must be chosen 
in such a way that 

U> (m (Cm{9) \{s}) =U g . (106) 

This induces the identification, 

Us({s}e IM ) = Us, g (107) 

of the byproduct operators. Now, the order of the desired 
unitary operation U g and the byproduct operator Us. g is 



as required in l|68[) . With adaptive measurement bases 
the effect of the randomness introduced by the measure- 
ments can be counteracted. What remains is a random 
byproduct operator which does not affect the determin- 
istic character of a QC c -computation and which is ac- 
counted for in the post-processing of the measurement 
results. 

In subsequent sections we will illustrate in a number 
of examples how Theorem ^ is used to demonstrate the 
functioning of quantum gate simulations on the QC C , and 
how the strategies for adapting the measurement bases 
are found. 



G. Function of CNOT-gate and general one-qubit 
rotations 

In this section, we demonstrate that the measurement 
patterns which we have introduced do indeed realize the 
desired quantum logic gates. 

The basis for all our considerations is the set 
of eigenvalue equations fulfilled by the cluster states. 
Therefore let us, before we turn to the realization of the 
gates in the universal set, describe how the eigenvalue 
equations can be manipulated. Equations are not the 
only eigenvalue equations satisfied by the cluster state. 
Instead, a vast number of other eigenvalue equations can 
be derived from them. 

The operators may for example be added, mul- 
tiplied by a scalar and multiplied with each other. In 
this way, a large number of eigenvalue equations can be 
generated from equations QJ. Note, however, that not 
all operators generated in this way are correlation oper- 
ators. Non-Hermitian operators can be generated which 
do not represent observables, yet will prove to be useful 
for the construction of new correlation operators. 

Furthermore, if quantum correlation operator K for 
state \4>) commutes with measured observable r"i ■ 
the correlation will still apply to the measured state. 
More specifically, if the state |</>) satisfies the eigenvalue 
equation K\(f>} = A|0) and [K, fi - a] = 0, then the 
state resulting from the measurement, Pj* |</i), where 
Psf — 1+( -~ 1 - > 2 ' T ''' cr — , satisfies the same eigenvalue equa- 
tion since A [P s ( f \4>)\ = [P^K\(j))} = K[P$ \ (f>)\. Thus the 
correlation K is inherited to the resultant state, P S j \(f>). 

To demonstrate and explain the measurement patterns 
realizing certain quantum gates, the program is as fol- 
lows. First, from the set of eigenvalue equations which 
define the cluster state \4>)c(g)> we derive a set of eigen- 
value equations which is compatible with the measure- 
ment pattern on Cm. Then, we use these to deduce the 
set of eigenvalue equations which define the state \tjj)c{g)j 
where the qubits in Cm have been measured. Thus we 
demonstrate that the assumptions for Theorem^ that 
is the set of equations (|82|) , are satisfied with the appro- 
priate unitary transformation U. Third, Us is obtained 
from equation Q84f) as a function of the measurement re- 
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suits. The order of U and Us is then interchanged and, 
in this way, the temporal ordering of the measurements 
becomes apparent. 

1. Identity gate 



the odd sub-lattice, we do not use it in the universal set 
of gates. Nevertheless, this realization of the Hadamard 
transformation can be a useful tool in gate construction. 
For example, we will use it in Section Hi G 41 to construct 
the realization of the z-rotations out of the realization of 
x-rotations. 



As a simple example, let us first consider a gate which 
realizes the identity operation 1 on a single logical qubit. 

For the identity gate Ci , Cm and Co each consist of a 
single qubit, so labeling the qubits 1, 2 and 3, 1 € Cj, 
2 e Cm and 3 e Co- The pattern M(t) corresponds to a 
measurement of qubit 2 in the a x basis. 

Let \<fi)c(i) be the cluster state on these three qubits. 
The state is defined by the following set of eigenvalue 
equations. 

4 X) 4 2) \4>)c ( i) =\4>)c { X), (108a) 
cfV^af \<t>) C (t> =\<f>)c m , (108b) 
4 2) 4 3) \4>)c(i) =\</>}c W - (108c) 

After the measurement of qubit 2, the resulting state 
of the cluster is 

IV-)c(i) = p^sMm, (109) 

where s 2 € {0, 1}, and = hthlTl^ , 

(2) (2) 

PxJ 2 an d °x obey the following relation, 

p£lA ] = {-^ S2 Pi% (no) 

Applying Pjfj 2 to both sides of equation (|108b() . and us- 
ing equation 1110(1 . one obtains for \ip}c{i)i defined in 
equation i|lU9|) . 

<#M 3) IV>>C(l) = (-l) S2 H')c(i)- (HI) 
Also from equations QlOSajl and IjlOScjl we have 

4 X) 4 3) |$ e(1) = \4>) C (i). (112) 

(2) 

Applying PxJ 2 to both sides of this equation gives 

<4 X) <4 3) IV>)e(i) = l^)c(i). (H3) 

Now, since qubits 1 and 3 represent the input and out- 
put qubits respectively, the assumption of Theorem ^ 
equation (|82|l . is satisfied for U = 1. The byproduct op- 
erator Us, is obtained from equation l|84[l . and we find 
that the full unitary operation realized by the gate is 
U = 1 < 2 erf 1 = <J S X 2 erf 1 1. 

Also note that a wire with length one (Ci(H) = 1, 
Cm{H) = 0, Co{H) = 2), i.e. half of the above elemen- 
tary wire, implements a Hadamard transformation. As 
in this construction the input- and output qubits lie on 
different sub-lattices of C, one on the even and one on 



2. Removing unnecessary measurements 

In larger measurement patterns, whenever pairs of ad- 
jacent a x - qubits in a wire are surrounded above and 
below by either vacant lattice sites or a z -measurements, 
they can be removed from the pattern without changing 
the logical operation of the gate. This is simple to show in 
the case of a linear cluster. Consider six qubits, labelled 
a to /, which are part of a longer line of qubits, prepared 
in a cluster state. Four of the eigenvalue equations which 
define the state are 

ai c Vi<M%) c = Me, 

Suppose, a measurement pattern M on these qubits 
contains measurements of the observable a x on qubits c 
and d. Measurements in the a x basis can be made before 
any other measurements in M. If these two measure- 
ments alone are carried out, the new state fulfills the fol- 
lowing eigenvalue equations, derived from equation l|114|) 
in the usual way, 

a^^a^)c = {-l) s ^)c, 

a^CT^cr^)c = {-l) s ^)c- 

The resulting state is therefore a cluster state from which 
qubits c and d have been removed, and b and e play the 
role of adjacent qubits. Thus, the two measurements have 
mapped a cluster state onto a cluster state and thus do 
not contribute to the logical operation realized by M 1 
which, in the case where both s c and Sd equal 0, is com- 
pletely equivalent to the reduced measurement pattern 
A4' , from which these adjacent a x measurements have 
been removed. 



3. One-qubit rotation around x-axis 

A one-qubit rotation through an angle a about the x- 
axis U x [a] = cxp[— ia/2a x ] is realized on the same three 
qubit layout as the identity gate. Labeling the qubits 1, 

2 and 3 as in the previous section, 1 = Cj, 2 = Cm and 

3 = Co- The measurement pattern M(U X ) consists of a 
measurement, on qubit 2, of the observable represented 
by the vector r xy {rf) = (cos(j]), sin (77), 0), 

r xy (r]) ■ v = cos r/a x + smi]cr y = U z [i]]a x U z [-r]}, (116) 



18 



whose eigenstates lie in the a?-y-plane of the Bloch sphere 
at an angle of rj to the x-axis. 

The cluster state \4>)c{u x ) IS defined by equa- 
tions l|l()8|) . After the measurement of A4(U X ), the re- 



sulting state is \ip) C {u x ) = P X y {v )\<t>)c{u x ) - xy(ri) 

2 2 r:cy ^' a . To generate an eigenvalue equation 
whose operator commutes with r xy (ri) ■ a we manipulate 
equation HU8cfl in the following way, 

<7 { J S) <r£ ) \<t>)c(u*) = \<t>)c(u x ) (H7) 
^%}c(^) = 4 3) l0>c(t/,) 
[4 2) -4 3) Mc (C /,) = 
.-. exphi(0/2)[af - a^McdJ*) = \<t>)c{U*) (118) 

where the last equation is true for all 9 € [0,27r]. This 
takes a more useful form, if we write it in terms of one- 
qubit rotations, 

uP[9]U^[-9M c{Uo:) = \^) c{Ux) . (119) 



where P 



(2) 



i.e. 
i.e. 



We use this and equation l|108b|) to construct the fol- 
lowing eigenvalue equation for \4>)c(Ux)i 



(120) 



Applying P xy ( n ),2 to both sides, we obtain the following 
eigenvalue equation for \ifi}c(Ux)> 
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In the same way as for the identity gate we also apply 
the projector to an eigenvalue equation generated from 
equations (|108a|) and H108c[) to obtain 



Wc(u»)=v£ ) v { * ) mc(u«) 



a^U^[- V }a^U^[ V }\^ c(Ua 



(122) 



and thus we see that equation 1821) is satisfied for U = 
U x {—rj) and /7s = cf 1 ^ 2 . Interchanging the order of 
these operators is not as trivial here as for the identity 
gate. When a z is propagated through U x (rj) the sign of 
the angle is reversed, so we find that the gate operation 
realized by this M(U X ) in the QC C is 



u g = u x [(-i) si (-v)} 



(123) 



The sign of the rotation realized by this gate is a function 
of si, the outcome of the measurement on qubit 1. This 
is an example of the temporal ordering of measurements 
in the QC C . In order to realize £4, [a] deterministically, 
the angle of the measurement, 77, on qubit 2 must be rj = 
(— l) Sl (—a), thus this measurement can only be realized 
after the measurement of qubit 1. 



I M(U g ) I = I M(HUgH) 



FIG. 5: Useful identity for the realization of the rotation U z [a] 
as the sequence H U x [a] H. 



4- Rotation around z-axis 

The measurement pattern for a rotation around the 
z-axis U z [/3] = exp[— if3/2a z ] is illustrated in Fig. [3 It 
requires 5 qubits for its realization. 

The measurement layout A4(U Z ) is similar to the ro- 
tation about the x-axis, except for two additional a x 
measurements on either side of the central qubit. The 
simplest way to understand this gate is regard it as the 
concatenation U z [a) = H U x [a] H . The Hadamard trans- 
formations may be realized as wires of length one, see 
Section Til G II Thus, the measurement pattern of the z- 
rotation is that of the x-rotation plus one cluster qubit on 
either side measured in the eigenbasis of a x , as displayed 
in Fig El 

The explanation in terms of eigenvalue equations 
obeyed by cluster states is as follows. Let us label 
the qubits 1 to 5. The cluster state \(j))c(u r ) is defined 
by eigenvalue equations of the usual form. If qubits 2 
and 4 are measured in the a x basis, the resulting state 

W)c(u z ) = Px%P x %\4>)c(u z ) fulfills the following set of 
eigenvalue equations 



c(u m ) = W)c<p.), (l 24a ) 
c(t/,) = (~1) S2 |0W)> (124b) 
^ 3 M% W) = (-IH^W)- (124c) 



This set of equations is analogous to equations (|108fl . 
except for the different eigenvalues and that the input 
and output qubits x- and z-bases have been exchanged. 
From here on the analysis of the measurement pattern 
runs parallel to the previous section. 

One finds A4(U Z ) realizes the operation U z ((3) if the 
basis of the measurement on qubit 3 is chosen to be 
the eigenbasis of (( — 1) S2 (—/?)) ■ c?, where r X y(rfj is de- 
fined in equation i|116|) . Qubit 2 must thus be measured 
prior to qubit 3. The byproduct operator for this gate is 



5. Arbitrary Rotation 

The arbitrary Euler rotation can be realized by com- 
bining the measurement patterns of rotations around x- 
and z-axes by overlaying input and output qubits of ad- 
jacent patterns, as described in section lll Dl This creates 
a measurement pattern of 7 qubits plus input and output 
qubits, labelled as in Fig. [SI with measurements of a x on 
qubits 3, 4, 6 and 7, and measurements in the x-y-planc 
at angles a, f3 and 7 on qubits 2, 5 and 8, respectively. 
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old labels 123456789 

EI ED □□ 



new labels 12 3 
If! removed - 



4 5 
Out 



FIG. 6: General rotation composed of two x-rotations and 
a z-rotation in between (Euler representation). In the QC C - 
realization pairs of adjacent cluster qubits measured in the 
CTz-eigenbasis may be removed from the measurement pattern. 



The unitary operation realized by these connected mea- 
surement patterns is, 

UxURotl^vX] =a>^C/ x [(-ir(-7)K 3+s ^+ S6 

(125) 

As we have shown above, adjacent pairs of a x measure- 
ments can be removed from the pattern without changing 
the operation realized by the gate. The operation real- 
ized by this reduced measurement pattern is obtained by 
setting the measurement results from the removed qubits 
to zero, S3, S4, sq, 37 = 0. After relabelling the remaining 
qubits in the measurement pattern 1 to 5, we obtain 



<7?o2»cu(-ir(-«)] 



(126) 



Propagating all byproduct operators to the left hand side 
we find the unitary operation realized by the measure- 
ment pattern is 



URodtvX] =C4H-i) sl+S3 7]£4H-i) S2 /3] 



J-S2+S4 0-S1+S3 



(127) 



One 



with byproduct operator Uy. 
finds that, to realize a specific rotation URot[t;,T), £] = 
^[Cl^W^nlfli * ne angles a, /?, 7 specifying the mea- 
surement bases of the qubits 2,3, and 4 are again de- 
pendent on the measurement results of other qubits. 
We see that a = (-l) Sl (-£), P = (-l) S2 (-r?), 7 = 
(— l) Sl+S3 (— £). To realize a specific rotation determin- 
istically, qubit 2 must thus be measured before qubits 3 
and 4, and qubit 3 before qubit 4, in the bases specified 
in Section HEEI 



6. Hadamard- and ir/2-phase gate 

The Hadamard- and the 7r/2-phase gate have the prop- 
erty that under conjugation with these gates Pauli oper- 
ators are mapped onto Pauli operators, 



Ha z H^ 



0~z 



(128) 



(129) 



and 

U M [ir/2]a x U,[ir/2]^ = o y 
U z [n/2}cT z U z [n/2}l = a z 

from which the propagation relations (|74|) follow. Re- 
lated to this property is the fact that these two special 
rotations may be realized via a x - and a y -measurements. 
Such measurement bases need not be adapted to previ- 
ously obtained measurement results and therefore, while 
these rotations might be realized in the same way as any 
other rotation, there is a more advantageous way to do 
so. 

To realize either of the gates we use again a cluster 
state of 5 qubits in a chain C(H). Let the labeling of the 
qubits be as in Fig. |3i and e, i.e. qubit 1 is the input- 
and qubit 5 the output qubit. 

A cluster state \4>)c{h) obeys the two eigenvalue equa- 
tions 

^(1)^(3)^(4)1 



)C{H) 



1C{H) 
)C(H) 



) c{ h) = K^K^K^\4>) C{H) 



(130) 



)C(H) 



When the qubits 2, 3 and 4 of this state are measured 
in the tr y -eigenbasis and thereby the measurement out- 
comes S2, S3, S4 € {0, 1} are obtained, the resulting state 
\'4>)c(H) obeys the eigenvalue equations 



(131) 



)C(H) = (-l) S3+Sl |0>C(ff), 
)C(H) = (-1) S2 + S3 |0>CW 

From equation (|128fl we see that the correlations 11311) 
are precisely those we need to explain the realization of 
the Hadamard gate. Using Theorem \I\ we find that by 
procedure|21with measurement of the operators a x \ Oy , 
<7y 3 ^ and tTy 4 ^ a Hadamard gate with a byproduct operator 
as given in H31f) is realized. 

A cluster state \4>)c(u z [n/2]) of a chain of 5 qubits obeys 
the eigenvalue equations 



)C(U z {tt/2]) 



A-(i)A-(3)a-Wa-(6)| 



p)c(u z [*m. 

}c(U z [n/2]) 



tf< 2 >#<%W [)r / 2]) 



(132) 



/C(U z [n/2]) 

)C(U,[n/2]] 

When the qubits 2, and 4 of this state are measured in the 
a x - and qubit 3 is measured in the ery-eigenbasis, with the 
measurement outcomes S2, S3, S4 <E {0, 1} obtained, the 
resulting state \^)c{u z W/2}) obeys the eigenvalue equa- 
tions 

.W„( 5 )i 



<7i CT y 



'V \v)c(U z [n/2]) = (-l) S3+S4+1 |'0}c{C/ z [7r/2]), 



(133) 

Using Theorem^ we find that by procedure f5] with mea- 



surement of the operators a x \ <j x 2 \ and cr'x 1 a n/2- 
phase gate is realized, where the byproduct operator is 
given by (|31|l . 



.(4) 
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1 2 3 4 5 6 7 

c - H □ H H H □ H c « 

^ □ □ □ □ H □ H l »- 

9 10 11 12 13 14 15 

FIG. 7: Pattern of correlation centers representing the eigen- 
value equation 1134at , 



7. The CNOT gate 

A measurement pattern which realizes a CNOT gate 
is illustrated in Fig. [21 Labeling the qubits as in Fig. [21 
we use the same analysis as above to show that this mea- 
surement pattern does indeed realize a CNOT gate in the 
QC C . 

Of the cluster C(CNOT) on which the gate is real- 
ized, qubits 1 and 9 belong to C/, qubits 7 and 15 belong 
to Co and the remaining qubits belong to Cm- Let \<j>) 
be a cluster state on C(CNOT), which obeys the set of 
eigenvalue equations 

From these basic eigenvalue equations there follow the 
equations 



= K {1) K {5) K {15) 

U X Uy Uy Uy Up Uy U ffl 

= K {2) iC (3) K {5) K {6) 



= CT (1) CT (2) (3) (5) (6) (7) 

U Z Uy Uy V y \J y U % 



lo) 



K (5) K (6) K (8) K (10) K (12) K (14) | 0) 

a (5) (6) (7)^(8)^(9)^(10) (12) (14) (15) 
w y w y w z y u z u x y x u z 



(134a) 
(134b) 
(134c) 
'}• 

(134d) 



Subsequently we will often use a graphic representation 
of eigenvalue equations like i|134afl - (|134d|l . Each of these 
equations is specified by the set of correlation centers q 
for which the basic correlation operators K^ q ' (J2J) enter 
the r.h.s. of the equation. While the information con- 
tent is the same, it is often more illustrative to display 
the pattern of correlation centers than to write down the 
corresponding cluster state eigenvalue equation. As an 
example, the pattern of correlation centers which repre- 
sents the eigenvalue equation i|134af> is given in Fig. [7| 

If the qubits 10, 11, 13 and 14 are measured in the 
a x - and the qubits 2, 3, 4, 5, 6, 8 and 12 are measured 
in the (7j,-eigenbasis, whereby the measurement results 
S2 - sq, S8: s io - s i4 are obtained, then the cluster state 
eigenvalue equations i|134a|) - (I134dll induce the following 



eigenvalue equations for the projected state l^) 
o£> <rjp 4 15) = (_i) 1 +«s+«*+*«+«»+«i3|^) j (135a) 



^ 9) 4 15) IV>) 



(_l)«a+«3+«5+«6|^ ( 135b ) 

(-1) S11+S13 |V), (135c) 

(_2)S 5 +S 6 +S S +Sio+Sl2+Sl4 1^/,) 

(135d) 



Therein, qubits 1 and 7 represent the input and out- 
put for the control qubit and qubits 9 and 15 represent 
the input and output for the target qubit. Writing the 
CNOT unitary operation on control and target qubits 
CNOT{c,t), we find 



CNOT{c, t)aV 


CNOT(c, t) 


u x u x i 


(136a) 


CNOT(c, t)af' 


CNOT(c,t) 


u z , 


(136b) 


CNOT(c, t)ag ] 


CNOT(c, t) 




(136c) 


CNOT(c, t)<j[ t] 


CNOT(c,t) 




(136d) 



Comparing these equations to the eigenvalue equations 
(|135a|) to 1135d|) . one sees that M. does indeed realize a 
CNOT gate. Furthermore, after reading off the operator 
Us using equations (|S2"|) and i|S4*)> and propagating the 
byproduct operators through to the output side of the 
CNOT gate, one finds the expressions for the byproduct 
operators, reported in equation l|24|) . 



H. Upper bounds on resource consumption 

Here we discuss the spatial, temporal and operational 
resources required for the QC C and compare with re- 
source requirements of a network quantum computer. 

To run a specific quantum algorithm, the QC C requires 
a cluster of a certain size. Therefore the QC c -spatial re- 
sources S are the number of cluster qubits in the required 
cluster state | </!))c, i.e. S = \C\. The computation is driven 
by one-qubit measurement only. Thus, a single one-qubit 
measurement is one unit of operational resources, and 
the QC C - operational resources O are defined as the total 
number of one-qubit measurements involved. The oper- 
ational resources O are always smaller or equal to the 
spatial resources S, 



0<S 7 



(137) 



since each cluster qubit is measured at most once. As for 
the temporal resources, the QC c -logical depth T is the 
minimum number of measurement rounds to which the 
measurements can be parallelized. 

Let us briefly recall the definition of these resources in 
the network model. The temporal resources are specified 
by the network logical depth T q \ n , which is the minimal 
number of steps to which quantum gates and readout 
measurements can be parallelized. The spatial resources 
5qi n count the number of logical qubits on which an al- 
gorithm runs. Finally, the operational resources O q i n are 
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the number of elementary operations required to carry 
out an algorithm, i.e. the number of gates and measure- 
ments. 

The construction kit for the simulation of quantum 
logic networks on the QC C shall contain a universal set 
of gates, in our case the CNOT gate between arbitrary 
qubits and the one qubit rotations. Already the next- 
neighbor CNOT with general rotations is universal since 
a general CNOT can be assembled of a next-neighbor 
CNOT and swap gates which can themselves be com- 
posed of next-neighbor CNOTs. However, in the follow- 
ing we would like to use for the general CNOT the less 
cumbersome construction described in Section llV CI For 
this gate, the distance between logical qubits, i.e. be- 
tween parallel qubit wires, is 4. The virtue of this gate is 
that it can always be realized on a vertical slice of width 6 
on the cluster, no matter how far control and target qubit 
are separated. A slice of width 6 means that the distance 
between an input qubit of the gate and the correspond- 
ing input of the consecutive gate is 6 lattice spacings. 
This general CNOT gate determines the spatial dimen- 
sions of a unit cell in the measurement patterns. The 
size of this unit cell is 4 x 6. The other elementary gates, 
the next-neighbor CNOT and the rotations arc smaller 
than a unit cell and therefore have to be stretched. This 
is easily accomplished. The next-neighbor CNOT as dis- 
played in Fig. [2^ has a size of 2 x 6 and is extended to 
size 4 x 6 by inserting two adjacent cluster qubits into 
the vertical bridge connecting the horizontal qubit lines. 
The general rotation as in Fig. |2t> has width 4 and is 
stretched to width 6 by inserting two cluster qubits just 
before the output. 

Concerning the temporal resources we first observe 
that we can realize the gates in the same temporal or- 
der as in the network model. To realize a general CNOT 
on the QC C takes one step of measurements, to realize 
a general rotation takes at most three. For the network 
model we do not assume that a general rotation has to 
be Euler-decomposcd. Rather we assume that in the net- 
work model a rotation can be realized in a single step. 
Thus the temporal resources of the QC C and in the net- 
work model are related via 

T < 3T q i n . (138) 

As for the spatial resources, let us consider a rectan- 
gular cluster of height h and width w on which the qubit 
wires are oriented horizontally, with the network regis- 
ter state propagating from left to right. As the logical 
qubits have distance 4, the height of the cluster has to 
be h = 4iS q i n — 3 where S q \ n is equal to the number n of 
logical qubits. Further, the number of gates in the cir- 
cuit is at most SqinTqin because, in the network model, 
in each step at most S q \ n gates can be realized. On each 
vertical slice of width 6 on the cluster there fits at least 
one gate such that -taking into account an extra slice 
of width 1 for the readout cluster qubits- for the width 
holds w < 6S q \ n T q \ u + 1. With S = hw one finds that 



In a similar way, a bound involving the network opera- 
tional resources can be obtained. The spatial overhead S 
and the operational overhead O per elementary network 
operation is < 2AS q \ n if this operation is a unitary gate 
from the universal set described before, and is equal to 
one if this operation is a readout measurement. Thus, we 
also have 



S < 24 q i n S'qln, 

O < 24 q i n 5 qln . 



(140) 



S < 24S' q i n 2 Tqi n . 



(139) 



The purpose of this section was to demonstrate that 
the scaling of spatial and temporal resources is at worst 
polynomial as compared to the network model. In Q 
it has been shown, as stated in Section IIII Al that the 
required classical processing increases the computation 
time only marginally (logarithmically in the number n of 
logical qubits) and thus there is no exponential overhead 
in either classical or quantum resources. 

The upper bounds in {T3g|) . ljT39"|) and ((TiUj) should 
not be taken for estimates. For algorithms of practical 
interest the required resources usually scale much more 
favorably and there do not even have to be overheads 
at all. This is illustrated for the temporal complexity 
of Clifford circuits in Section ITTT1 and in the examples of 
Section Hvl A spatial overhead always exists. However, 
this is compensated by the fact that the operational effort 
to create a cluster state is independent of the cluster size. 



I. Quantum circuits in the Clifford group can be 
realized in a single step 

The measurement bases to realize the Hadamard- and 
the 7r/2-phasc gate need not be adapted since only op- 
erators <j x and tjy are measured. The same holds for 
the realization of the CNOT gate, see Fig. |21 Thus, all 
the Hadamard-, 7r/2-phase- and CNOT-gates of a quan- 
tum circuit can be realized simultaneously in the first 
measurement round, regardless of their location in the 
network. In particular, quantum circuits which consist 
only of such gates, i.e. circuits in the Clifford group, can 
be realized in a single time step. As an example, many 
circuits for coding and decoding are in the Clifford group. 

The fact that quantum circuits in the Clifford group 
can be realized in a single time step has previously not 
been known for networks. The best upper bound on the 
logical depth that was known previously scales logarith- 
mically with the number of logical qubits 01 ■ 

Note that, as stated by the Gottesman-Knill-Thcorcm 
[l9j , there is no need for fast Clifford circuits if the quan- 
tum output is measured in a Pauli basis because these 
circuits can be simulated efficiently classically. However, 
the purpose of this section is to point out that the whole 
Clifford part of any quantum circuit can be performed in 
a single time step. We will discuss this point further in 
Section IITTbI 

Here we find a first aspect of QC c -computation which 
is not adequately described within the network model, 
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and with this observation we conclude the discussion of 
the QC C simulator of quantum logic networks. 



III. COMPUTATIONAL MODEL UNDERLYING 
THE QC C 

A. Processing of information 

In the network model of quantum computation one 
usually regards a quantum register as the carrier of in- 
formation. The quantum register is prepared in some in- 
put state and processed to some output state by applying 
a suitable unitary transformation composed of quantum 
gates. Finally, the output state of the quantum register 
is measured by which the classical readout is obtained. 

For the QC C the notions of "quantum input" and 
"quantum output" have no genuine meaning if we re- 
strict ourselves to the situation where the input state 
is known. As stated before, Shor's factoring algorithm 
[l5| and Grover's search algorithm are both exam- 
ples of such a situation. In these cases the final result 
of any computation -including quantum computations- 
is a classical number. In a QC C -computation this num- 
ber is extracted from the outcomes of all the one-qubit 
measurements. The entire computation amounts to just 
measurements of the cluster qubits in a certain order and 
basis. 

We have divided the set C of cluster qubits into sub- 
sets /, M and O to describe the QC C in terms of the 
network model. Such a terminology is not required for 
the QC C a priori. It is true that when a quantum logic 
network is realized on a cluster state there is a subset 
of cluster qubits which play the role of the output reg- 
ister. However, these qubits are not the final ones to 
be measured, but among the first (!). The measurement 
outcomes from all the cluster qubits contribute to the 
result of the computation. The qubits of O C C simu- 
late the output state of the quantum register and thus 
contribute obviously to the computational result. The 
cluster qubits in the set I <Z C simulate the fiducial in- 
put state of the quantum register and their measurement 
contributes via the accumulated byproduct operator on 
O. Finally, the qubits in the section M C C of the clus- 
ter whose measurements simulate the quantum gates also 
contribute via the byproduct operator. 

Naturally there arises the question whether there is any 
difference in the way how measurements of cluster qubits 
in /, O or M contribute to the final result of the compu- 
tation. As shown in Q, it turns out that there is none. 
This is why we can abandon the notions of quantum in- 
put, quantum output and quantum register altogether 
from the description of the QC C . Furthermore, quantum 
gates are not constitutive elements of the QC C ; these are 
instead one-qubit measurements performed in a certain 
temporal order and in a spatial pattern of adaptive mea- 
surement bases. In fact, the most efficient temporal order 
of the measurements does not follow from the temporal 
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FIG. 8: General scheme of the quantum computer via one- 
qubit measurements. The sets Qt of lattice qubits are mea- 
sured one after the other. The results of earlier measurements 
determine the measurement bases of later ones. All classical 
information from the measurement results needed to steer the 
QC C is contained in the information flow vector I(t). After 
the last measurement round tmax, I(imax) contains the result 
of the computation. 



order of the simulated gates in the network model. 

The general view of a QC c -computation is as follows. 
The cluster C is divided into disjoint subsets Qt C C with 
< t < t max , i.e. |Jt=p Qt = C and Qs H Qt = for all 
s =/= t. The cluster qubits within each set Q t can be mea- 
sured simultaneously and the sets are measured one after 
another. The set Qo consists of all those qubits for which 
no measurement bases have to be adjusted, i.e. those of 
which the operator a x , a y or a z is measured. In the sub- 
sequent measurement rounds only operators of the form 
cos tpcTx ± sin^cTj, are measured where \tp\ < 7r/2, (p ^= 0. 
The measurement bases are adaptive in these rounds, 
i.e. they are adapted to measurement results obtained in 
previous rounds. The measurement outcomes from the 
qubits in Qo determine the measurement bases for the 
qubits in Qi, which are measured in the second round, 
those from Qo and Qi together determine the bases for 
the measurements of the qubits in Q2 which are mea- 
sured in the third round, and so on. Finally, the result 
of the computation is calculated from the measurement 
outcomes obtained in all the measurement rounds. 

Now there arises the question of how complex the re- 
quired classical processing is. In principle it could be that 
all the obtained measurement results had to be stored 
separately and the functions to compute the measure- 
ment bases were so complicated that one would gain no 
advantage over the classical algorithm for the considered 
problem. This is not at all the case. If the network algo- 
rithm runs on n qubits then the classical data that the 
QC C has to keep track of is entirely contained in a 2n- 
componcnt binary valued vector,which we have called 
the information flow vector I(t) Q- The update of I(t) 
is a classical computation that is needed to adapt the 
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measurement bases of cluster qubits according to previ- 
ous measurement outcomes. These updates and the final 
identification of the computational result from I(i max ) 
are all elementary. 

Concerning the resources for the classical processing 
of the measurement outcomes in a QC c -computation, we 
point out that this processing increases the total time of 
computation only marginally 0. 

In summary, the formal description of the QC C is based 
on primitive quantities of which the most important ones 
are the sets Qt C C of cluster qubits defining the temporal 
ordering of measurements on the cluster state, and the 
binary valued information flow vector 1(f) which is the 
carrier of the algorithmic information. The reader who 
is interested in how this computational model arises and 
in its detailed description is referred to • 



B. Algorithms and graphs 

In this section we relate QC c -algorithms to graphs. We 
do this by considering non-universal graph states suited 
for the specific algorithm in question. For the QC C , the 
Clifford part of each algorithm can be removed. A math- 
ematical graph comprises all the information that needs 
to be kept the Clifford part. 

While the network formulation of a quantum algorithm 
is given as a sequence of quantum gates applied to a fidu- 
cial input state, the QC c -version of a quantum algorithm 
is specified by a measurement pattern on the universal 
cluster state plus the structure for the processing of 
the measurement outcomes. 

The measurement pattern is, in the simplest case, just 
a copy of the network layout to the substrate cluster 
state, imprinted by the measurements. As such it con- 
tains information about the precise location of the gate 
simulations and about the way the "wires" connecting 
the gates are bent around. These are all details of the 
realization of an algorithm but do not belong to the de- 
scription of the algorithm itself. Thus, the measurement 
pattern introduces a large amount of redundancy into the 
description of a QC c -algorithm. This redundancy may be 
reduced to a large extent by allowing for non-universal, 
algorithm-specific quantum resources. 

Clearly, at this point one has to specify how special the 
algorithm-specific resource is allowed to be. Obviously 
it would make no sense to take the quantum output of 
the entire network as the required quantum resource and 
to regard the subsequent readout measurements as the 
algorithm. Here, we allow for any graph state 0], 12 1|) 
as the quantum resource. Graph states are both easy to 
create and to describe. Every algorithm may be run with 
a graph state as the quantum resource since the cluster 
state is a particular graph state. 

To allow for an algorithm-specific graph state as the 
quantum resource of a QCc-computation reduces the re- 
dundancy of both the description and the realization of 
a quantum algorithm. This can easily be seen from the 



material presented in Section Hi CI All the cluster qubits 
q G C\C/v can be get rid of either by measuring them 
in the a z -eigenbasis or cquivalently by not placing them 
initially into their positions at all. The remaining state 
on the sub-cluster Cat is again a cluster state. Hence it 
is also a graph state. It is less redundant and no longer 
universal. 

But we can go further. Not only the qubits measured 
in the oveigenbasis may be removed from the cluster 
but instead all those qubits of which one of the Pauli 
operators a x , a y or cr z is measured, i.e. all the qubits 
which form the set Qq. The state of the unmeasured 
qubits that emerges after the measurement of the cluster 
qubits in Qo is again (local equivalent to) a graph state. 

This may be seen as follows. First note that the op- 
erators ai^ (£) beV (v'z'j which appear in (|2"T|) form 
a stabilizer of the state \4>{k])g- The generator of the 
stabilizer contains \C\ elements for a state of \C\ qubits. 
After all the qubits q £ Qq have been measured, the re- 
sulting state |^ , )c\Q °f the |C\Qo| unmeasured qubits is 
again described by a stabilizer of the form 

\C\Qo\ x z 

(8) ""(4°) i*) C \ Qo = ±i*) CVQo (141) 

Va = l..\C\Q \, 

with two |C\Qo| x |C\Qo|-matrixes X and Z, for which 
X a .i, Z a ^ G {0, 1}. The \C\Qo\ x 2|C\Q |-compound ma- 
trix (X\Z) is called the generator matrix of the stabilizer 
for \ty)c\Q - The state | v I / )c\Qo 1S uniquely determined 
by the generator of its stabilizer. 

The state | v E , )c\q can thus be regarded as a 
[|C\Qo|j 0, d]-stabilizer code, with the distance d not spec- 
ified. Whether a code with only one code word that en- 
codes qubits should be regarded as a code in the sense 
of coding shall not concern us here. For the present pur- 
pose, it is important to note that the state \^)e\Q fulfills 
the assumptions of Theorem 1 in [2(j ■ The cited theorem 
states that any stabilizer code over the alphabet A = ¥ p m 
is [local unitary] equivalent to a graph code. If we special- 
ize to our case, A = F 2 2, we find that the state \^)c\q 
specified in (|141f> is local unitary equivalent to a graph 
state \(/>{K}} G{ c\Q ,E C y Qo ) 63- 

That is, the state \^)c\q needed for the computation 
can be obtained from a graph state \<I>{k}) g(c\q ,e C \ Q ) 
via local unitary transformations. Subsequently in the 
process of computation, the qubits of \^)c\Q are mea- 
sured in their given temporal order and in the appro- 
priate adapted bases. An alternative way to proceed is 
to use the graph state \H k })g(c\Qo,e CKQo ) directly, only 
modifying the measurement bases instead of performing 
the local rotations prior to the measurements. Thus, in a 
QC c -computation with a special graph state as the quan- 
tum resource and the first measurement round omitted, 
the way of processing the classical information is the same 
as in a QC c -computation with a universal resource and 
the first measurement round performed. 
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The graphs associated with states (1141(1 are in general 
not unique [2(| ■ A constructive way to obtain graphs on 
C\Qo from G(C,Ec) and the measurement bases of the 
qubits in Q has been described in |2lj . 

Now note that the measurement of the qubits in Qq 
realize the Clifford part of a quantum circuit. The fact 
that we can reduce the quantum resource by these qubits 
means that we can remove from each quantum algorithm 
its Clifford part. This represents, in a way, an exten- 
sion to the Knill-Gottesman-Theorem stating that 
a quantum computation that consist only of quantum 
input state preparation in the computational basis, uni- 
tary gates in the Clifford group, measurement of obscrv- 
ables in the Pauli group, and gates in the Clifford group 
conditioned on the outcomes of such measurements, may 
be simulated efficiently classically and thus requires no 
quantum resources at all. 

With only a single non-Clifford operation in the cir- 
cuit, such as a one-qubit rotation about most axes and 
angles, the efficient classical formalism upon which the 
Gottesman-Knill theorem rests can no longer be applied. 
The QC c -construction, on the other hand, is not affected 
by this. Each quantum network algorithm in question 
may be reduced by its Clifford part. Only the non- 
Clifford gates require quantum resources. The price is 
that the universal quantum resource, the cluster state, is 
changed into a non-universal, algorithm-specific resource 
-a graph state l|21(l- on fewer qubits. The Clifford part of 
the network algorithm specifies the corresponding graph. 

In conclusion, instead of describing a quantum algo- 
rithm as a network of gates applied to some fiducial input 
state, a quantum algorithm may (arguably more effec- 
tively) be characterized by a graph specifying the quan- 
tum resource and the structure Q for the processing of 
the measurement outcomes. 



IV. EXAMPLES OF PRACTICAL INTEREST 

A. Multi-qubit swap gate 

A multi-qubit swap gate is an n-qubit generalization of 
the two-qubit swap gate. It reverses the order of the n- 
qubits, interchanging qubit i with n+ 1 — i, i = 1, 2, .., N. 
This can be realized in a simple way on the QC C , as 
shown in Fig. [5Jt. The measurement pattern M. on Cm 
consists of a square of a x measurements, with sides of 
In — 1 cluster qubits. The input qubits are, simultane- 
ously with the qubits in Cm, also measured in the o~ x - 
eigenbasis. 

It can be verified using the methods introduced above 
that realizing M. leads to correlations between the ith 
input qubit and the n + 1 — ith output qubit. Here, we 
discuss the four-qubit swap as a particular example. 

After the <r x -measurements of the qubits in Cm wc 
obtain for the projected state IV')c(swap) the eigenvalue 
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FIG. 9: The multi-qubit swap gate, a) measurement pat- 
tern to realize the swap gate, b) Correlation centers for 
two correlations of the projected state |'i/'}c(swap) as inherited 
from correlations of </>)c(swap) ■ The correlation ai'' 1 ^ <ri ' 4 ' of 
| / i/>}c(swap) stems from the product correlation for |</>)c(swap) 
with the centers a of basic correlation operators de- 
noted by "o". The centers of the initial correlation, which 
after the measurements induces the correlation ai 1 ' 3 ^ o"z 2 ' of 
l^)c(swap), are denoted by 



equations 








hA)c(swap) = 


(-l) Ax ' 1 |V')c(swap) 




l^)c(swap) = 


(-l) Ax ' 2 |'0)c(swap) 


a ^3) (0,2) 

^ X W X 


IV')c(swap) = 


(-l) Ax ' 3 |'0)c(swap) 


^ X X 


|V'}c(swap) = 


(-1) X ' 4 |'0)c(swap) 


^ /!l) 4°' 4) 


hA)c(swap) = 


(-1) Inswap) 




V')c(swap) = 


(-1) "' 2 |V , >C( S wap) 




h/')c(swap) = 


(-l)^ 8 Inswap) 


CT (/,4) CT (0,1) 

w z u z 


V')c(swap) = 


(-1) zA \lp) C (sw a p) 



(142) 



Therein, the parameters \k.x,^k,z €E {0,1} depend lin- 
early on the measurement outcomes {s(j,j)}- Therein, 
i is the value of the x- and j the value of the y- 
coordinate of the respective qubit site. For example, 

-W = «(1,2) + s (2,3) + s (3,4) + S(4,5) + s (5,6) + s (6,7) m °d 2. 

The eigenvalue equations ll 142(1 can be derived from 
corresponding eigenvalue equations for the cluster state 
I0)c(swap) on the cluster C(swap). The required initial 
correlations are products of the basic correlation oper- 
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ators The way to obtain the equations l|142l) is 

rather straightforward and therefore we omit the detailed 
derivations. In Fig. Eh) two examples for the composition 
of product correlation operators from basic correlation 
operators are illustrated. The first line of l|142|) . 

vVVoVA Inswap) - Inswap), 

for example, is derived from the eigenvalue equation 

Ifewap) =^ (C ^ l) |0)c(swap), (143) 

with 

a£C x ,i 

and C x ,i = {(1,1), (1,2), (2, 3), (3, 4), (4, 5), (5, 6), (6,7), 
(0,4)}. Multiplying it out we find that all operators a z 
cancel and that 

= ' a x a) . (145) 

It is now easy to see that after the a x measurements 
of the qubits in Cm there remains a strict ai <?x ' - 
correlation for the state | V)c ( swap ). A similar construc- 
tion can be given to obtain the ai 1 ' 1 ^ eri°' 4 ^-correlation. 

With the eigenvalue equations l|142fl the assumptions 
of thcorcmnai'c fulfilled and thus via the described mea- 
surement pattern a unitary operation U = SWAP is re- 
alized modulo a byproduct operator as specified in (|84fl . 
To exchange the order of the swap-gate £4wap and the 
byproduct operator Us the byproduct operator is con- 
jugated under J7 swap , as usual for gates in the Clifford 
group. 

B. Simulating multi-qubit Hamiltonians 

Here we display a gate which simulates the unitary 
evolution with U = cxp(— iH^t) of the quantum input 
for the multi-particle Hamiltonian 

H A =gaVaf)*f)a<f> (146) 

and arbitrary times t. In addition, the gate performs a 
swap which can be corrected for by a subsequent swap 
gate as described in Section TlV Al 

The procedure to realize the measurement pattern M. 
for Hamiltonian simulation, as shown in Fig. 1101 requires 
two rounds of measurements. In the first round all the 
<j x -measurements are performed. In the second measure- 
ment round, of the qubit (3,4) the operator 

r(3,4) • o = U z [(-l) x "'2ip] a x C/][(-l) A "2^] (147) 

is measured, where U z [a] = exp(— iota z / 2). Therein, the 
angle ip is given by 

ip = gt, (148) 



and Am & {0, 1}, which depends linearly on outcomes of 
measurements in the first round, will be specified below. 

To understand the functioning of the Hamiltonian sim- 
ulator let us first discuss the state \ip') on the cluster 
C(sim) after the first round of measurements. By argu- 
ments analogous to those used in Section II V Al the state 
\ip') obeys the following eigenvalue equations: 
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Further, the state obeys the eigenvalue equation 

CT (3,4) ^(0,1)^(0,2)^(0,3)^(0,4) |^/ ) = (-l)*^), (1 50 ) 

with A £ {0, 1} linear in the measurement outcomes of 
the first round. Equation i|15U|) can be easily verified with 
the pattern of correlation centers displayed in Fig. 110b . 
From (|15Uf> it follows that 

cxp (i6a^) U4(-l) x 9] |V0 = W) (151) 

for arbitrary angles 9, with 

Ui[a] = exp (-zaaf^a^^a^^af^) . (152) 

Equation l|151|) is now inserted in both the l.h.s. and 
r.h.s. of the equations (|149|) . For example, with the first 
equation from i|149|) one obtains 

= (U z [26]a x Ut\2e]f A) o^ 

[u i [-{-ife]a x ^ul[-(-ife]) ( °\^). 

(153) 

In the second measurement round the qubit (3, 4) is the 
only one left to be measured. As can be seen from (|153fl . 
if of the operator U z [28]a x Ul[26] of qubit (3,4) is mea- 
sured then the state \tp), into which the cluster qubits are 
projected after the second measurement round, obeys the 
eigenvalue equation 

(-1) A »,1+ S (3,4) = 

<r™ (c/ 4 [-(-l)^]4 4 ^l[-(-l)^]) (0) |# 

(154) 

If we carry out this procedure for all equations in l|149|) 
we find that the state that emerges after the second 
measurement round obeys the eigenvalue equations 

(u^^U^Ul)^ |V) = (-1) A -+ S <^> |V), 
a™ {u^ ap a^U avJ 4ul) (0) |V) = (-1) A " |V), 

(155) 
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for i = 1 . . 4 and with U4 written in short for 

c/ 4 [-(-i) A 0]. 

With the set of equations 1)155(1 the assumptions 18211 
of theorem n are fulfilled. With theorem ^ it follows that 
the measurement pattern displayed in Fig. 1101 realizes a 
unitary transformation 



U sim = U4-(-l) x 6]U swlip Uz, 
where the byproduct operator is given by 

A 

s (/,i)+ A x,i+S(3,4) 



(4<>) 



A* 



(156) 



(157) 



Finally, the order of the operators has to be exchanged. 
Note that C/ swap and U4 commute. From (|156|) one finds 



u^u swap u4-(-i) x+ zt=i**,i e], 



with 



(158) 
(159) 



Thus, in order to realize U^ip] with <p specified in 11481) 
we must choose 



(160) 



That is, in the second measurement round we measure 
on the qubit (3,4) the operator given in l|147|) . where 
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FIG. 10: Simulation of the Hamiltonian Hi as specified in 
eq. 11461 . a) measurement pattern, b) Correlation centers 
for additional correlation. Shaded squares (in b)) represent 
cluster qubits measured in adaptive bases. 



Am = 1 + A + A * 



mod 2. 



(161) 



The {\ x ,i}, {^z,i} and A depend linearly on the measure- 
ment outcomes {s(j.j)} obtained in the first measurement 
round. 

The sub-circuit we have described in this section simu- 
lates the unitary evolution according to a particular four- 
particle Hamiltonian in a two-step process of measure- 
ments. The time for which the simulated Hamiltonian 
acts is encoded in the basis of the measurement in the 
second round. 

The generalization of the simulation of the 4-particlc 
Hamiltonian H±, shown in Fig. 1101 to an arbitrary num- 
ber n of qubits, i.e. the simulation of the Hamiltonian 
H n = ®" = i &z\ is straightforward. 



C. CNOT between non-neighbouring qubits 
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The CNOT gate described in Section III G 71 operates 
on two logical qubits whose input qubits are adjacent to 
each other on the cluster. However, for universal quan- 
tum computation, one must be able to realize a CNOT 
gate between any two logical qubits. While this could be 
achieved using a combination of the CNOT gate, intro- 
duced above, and the swap gate, the width of the mea- 
surement pattern needed to realize this would grow lin- 
early with the separation of the two logical qubits. There 



FIG. 11: Measurement pattern for a CNOT gate between two 
logical qubits whose input and output qubits are not neigh- 
bors. Squares in light gray denote cluster qubits measured in 
the eigenbasis of a x , in dark gray of a y . Pattern (a) is for the 
case where the two qubits are separated by an odd number 
of logical qubits. Pattern (b) is for an even numbered sep- 
aration. The patterns can be adapted to any separation by 
repeating the section enclosed by the dashed line. The width 
of the pattern remains the same for all separations. 
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bin 
□ 



a. ni 



□ a out 



□ 

bom 




FIG. 12: This measurement pattern is one of the key compo- 
nents of the measurement pattern in Fig. 1111 It performs a 
conditional 7r-phase-gate and a swap-gate. 




FIG. 13: The measurement pattern in Fig. 1111 realizes the 
quantum logic circuit on the left hand side of this figure. This 
network is equivalent to the one on the right hand side, where 
the only gate realized is the CNOT between the two desired 
non- adjacent qubits. 



is, however, an alternative measurement pattern, which, 
at the cost of doubling the spacing between the input 
qubits on the cluster, has a fixed width. The measure- 
ment pattern is illustrated in Fig. llll for qubits separated 
by an odd and even number of logical qubits, respectively. 

This layout can be understood within the quantum 
logic network model. The "wires" for the logical qubits in 
between the target- and the control qubit are crossed us- 
ing the measurement sub-pattern illustrated in Fig. 112b . 
However, as well as swapping the qubits, this pattern 
also realizes the a controlled 7r-phase gate, also known as 
a controlled a z gate, illustrated in Fig. fT5b . 

The quantum logic circuit realized by the whole mea- 
surement pattern, illustrated on the left-hand side of 
Fig. EH uses these sub-patterns to swap the positions of 
adjacent qubits. This brings non-neighboring qubits to- 
gether so that a CNOT operation may be performed on 
them. 

The networks on the left and on the right of Fig. 1 131 act 
identically, and thus the measurement pattern displayed 
in Fig. realizes a distant CNOT-gate. 



D. Controlled Phase Gate 

Here, we give an example of another two-qubit 
gate which can be realized without decomposing it 
into CNOTs and rotations, the controlled phase gate 



1 2 3 4 x 
□ 

\-(lb) □ □ (O.b) 



In Out 

FIG. 14: Controlled phase gate with additional swap. 

£^CPg($)- This gate realizes the unitary operation 

Ucpg[9] = l (ab) + (e lf> - 1) |H)a6(ll|, (162) 

applied to the two qubits a and b. 

We can write this in terms of the following one- and 
two-qubit rotations, 

Ucpg[0] =e^uit ) [-e/2}U^[6/2]U^[e/2}, (163) 
where the two-qubit rotation is 

U^[6]=e^(-\e/2a^af)). (164) 

This representation is particularly convenient for find- 
ing the measurement pattern that realizes the gate, since 
rotations U z [6/2] and U zz [—8/2] are realized on the 
QC C in a simple natural way. The measurement pat- 
tern is illustrated in Fig. O in which the labelling of the 
qubits is also defined. 

We follow the same method as above, beginning with 
the eigenvalue equations of the cluster state \<j>)c on the 
qubits shown. The er^-mcasuremcnts can be considered 
first, using the methods already illustrated in this paper. 
The resultant state of the remaining qubits \ip'}, after 
this sub-set of the measurements has been carried out, is 
defined by the following set of eigenvalue equations. 

oMoM>oM>o*>°>W) = W), (165a) 
oMa^a^M^m = \n (165b) 

(165c) 

a {b,I) a {a,0)^ = (_ 1 ). (I ,s)+«( a , a) +«cs. 1 )|^') j 

(165d) 



and 



^(2,1)^(0,0)1^ = (-IJ-tM)^), ( 166a ) 

ai 2 'V4 b >°^') = (-iy™\i>'), (166b) 

CT (l,2) CT (a,0) CT (b,0)|^ = (_ 1 )«cs, I )+«Ca. a )+»(s,s)|^. 

(166c) 



As in section Hi G 31 eigenvalue equations are now gen- 
erated which commute with the remaining measurements 
in Jvl, namely the measurements of axy(cti) on qubits 
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i G {(2, 1), (1, 2), (2, 3)}. First, wc manipulate the equa- 
tions l|166[) such that, for example, the eigenvalue equa- 
tion l|166cjl attains the form 

Ui 1 ' 2) [£,]U ( z i°' a) ' i °' b)) [-(-1) S (3,1)+ S (2,2)+S(3,3),C] \tP') 

= W)- 

(167) 

Similar equations containing one-qubit rotations on 
qubits (2, 1) and (O, a), and (2, 3) and (O, b) are derived 
from the other equations of 1166|) in the same way. These 
equations are inserted into both sides of the eigenvalue 
equations l|165|) so that, using the method introduced 
above, we obtain a set of four eigenvalue equations for 
\ip') which induce a set of four eigenvalue equations for 



Using 1)1 70J) one finds the following result: To re- 
alize the controlled phase gate 1162|) together with a 
swap-gate, the observables (|168|l measured in the sec- 
ond round have to be chosen with the angles a^.i) ~ 

1+0(3,2) +«(1,3) 0/2, a (1 , 2 ) = (-l) S <l.l)+ S (2.2)+«(1.3)6)/2 

and a (2 , 3 ) = (-l) s (u)+ s (2,2) +1 fl/2. This realizes the 
gate U-eUcpg[8], where the byproduct operator Us gen- 
erated by the measurements may be read off from equa- 
tion fl7U|> . 

E. Quantum Fourier transformation 

To realize the quantum Fourier transform we simulate 
the quantum logic network given in Fig. 115b . The ar- 
rangement of the gates in this network is taken from 
[22^ . Note that in [23 it was demonstrated that the setup 
to perform a quantum Fourier transformation simplifies 
considerably in a situation where the output state is mea- 
sured right after the transformation. Here, however, the 
quantum Fourier transformation may constitute part of a 
larger quantum circuit and we do not measure its output 
state. 

As can be seen from Fig. 1151 the quantum Fourier 
transform consists of Hadamard gates and combined 
gates which perform a conditional phase shift and a swap. 
These gates have been discussed in Sections III Bl and 
IIVDI All that remains to do is put the measurement pat- 



the state \tp) that one obtains after the remaining mea- 
surements have been carried out. 

Specifically, in the second measurement round the 
qubits (1,2), (2,1) and (2,3) are measured. Of these 
qubits one measures the observables 

r a -a^ = {U z [a a ]a x U z [a^) {a \ (168) 

for a G {(1, 2), (2, 1), (2, 3)} and the {a a } specified below. 

The induced eigenvalue equations for the state are 
of the form of equation (|82() , and the unitary operation 
realized by the gate can be read off from them using 
theorem ^ The full unitary operation realized by the 
measurement pattern is 



(169) 



(170) 




FIG. 15: Quantum Fourier Transformation, a) Network for 
quantum Fourier transformation on four qubits, taken from 
I22jl . b) Component of the network shown in a) which per- 
forms a conditional phase- and a swap-gate. Specifically, the 
gate shown is Ucpg [27r/2 m ], i.e. U m = |0> <0| + e^ 2 "" |1)(1|. 



terns simulating these gates together, using the network- 
like composition principle described in Section III Dl 

In this way we obtain a measurement pattern in which 
there are adjacent cluster qubits in "wires" that are mea- 
sured in the a x -eigenbasis. As described in Section lll G 21 
such pairs of cluster qubits may be removed from the 
measurement pattern. Note, that by removing adja- 



u'Uv = u£V[-(-iy^ +s ^ + ^ 

/ ( a )~j S (l,l)+ S (2,2)+S(3,3) / ( 6 ^S(1,3)+S(2,2)+S(3,1) / ( Q ) \0 (/,„) +0 ( j >2 ) +S (2 3) , ^ s S(j, S ) +S(2,1) +0(1,2) 

such that after the order of the gate and the byproduct operator is reversed, U'U S = UsU, one obtains 

tt tt ( ( a )\ S(1 ' 3 > +;S < 2 ' 2 > +S ( 3 ' 1 > / (b)\ S (l.l)+ S (2,2)+S(3,3) / ( a )\ s (2,l)+ s (l,2)+S(I,b) / 0(/ )a ) +0(1,2) +0(2,3) 

UyU = \(Tx \ fi' kz <7z 



u£V[-(-iy^ +s ^ +s ^^ 
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In Out 

FIG. 16: QC c -realization of a quantum Fourier transforma- 
tion on four qubits. The cluster qubits displayed as framed 
squares are measured in adapted bases. For the labels see 
text. 



cent pairs of cr^-measured cluster qubits we have moved 
the Cy-measurements of the Hadamard transformations 
"into" the subsequent conditional phase gates, i.e. we re- 
moved a cluster qubit which was not from a wire. It can 
be easily verified that this is an allowed extension of the 
method described in Section Til G 21 Finally, one obtains 
the QC c -circuit displayed in Fig. ^] 

In this circuit, as in all the others, the adaptive mea- 
surements are of observables 
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FIG. 17: A measurement layout to realize a Toffoli phase gate 
with phase <f>. The qubits marked by black boxes are simul- 
taneously measured in adapted bases depending on previous 
measurement outcomes. 
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FIG. 18: Cluster state quantum correlations for the realiza- 
tion of [4>/Q, used in the Toffoli phase gate. 



U z [±r,]a x U z [±itf 



(171) 



with rj = 7r/4 for cluster qubits marked with "2" in 
Fig. ^3 rj = 7r/8 for qubits marked with "3" and r\ = 
ir/16 for the qubits marked with "4". The sign factors 
of the angles in i|171|) depend on the results of previous 
measurements . 

The QC c -circuit, shown in Fig. ^| for the case of four 
qubits, is straightforwardly generalized to an arbitrary 
number n of logical qubits. The temporal spatial and 
operational resources T, S and O are, to leading order 



T = n, S,0 = 2n 2 



(172) 



The corresponding network resources are T q i n = 2n, 
Sq\ n = n and O q i n = n 2 /2. Thus, the scaling of the 
QC C spatial resources is worse than in the network 
model, but the temporal and operational resources scale 
in the same way as the corresponding resources for the 
network. The QC c -simulation of the network displayed 
in Fig. 1 151 requires half as many time steps and four times 
as many operations, albeit only one-qubit operations. 



F. Multi-qubit controlled gates 

In this section we describe the realization of the Toffoli 
phase gate and the three-qubit controlled gate CARRY 
which wc will both need for the construction of the QC C - 
addcr circuit described in Section HV Gl 

The Toffoli phase gate is a three-qubit generalization 
of the two-qubit controlled phase gate. If all three qubits 
arc in the state |1), the state gains a phase of exp(i</>), 



while all other logical basis states remain unchanged by 
the gate, 

4offo C H t) M = l (ci ' C2,4) + (e # -l)|lll)c 1 ,c 3 , t (lll|. (173) 

Like the controlled phase gate it can be represented as a 
product of multi-qubit rotations, 



rr(ci,c 2 ,t) r,i _ rr( Cl ,c 2 ,t) 
^Toffoli m — u zzz 





7T(ci,C2) 


4>~ 


rr(ci.*) 


f 


_4_ 


^ 22 


4_ 


^ zz 


4_ 



[/(d) 



4 

(174) 



where we have dropped the global phase, and 

Utz c2,t) [a] = exp (-fo/2ai Cl) ( ri C2) ai t) ) is a three qubit 

generalized rotation. The two-qubit rotations U zz are as 
defined in Q164p . 

The way to convert the sequence ll 174(1 of generalized 
rotations into a measurement pattern is as in the ex- 
amples before. The measurement layout for the Toffoli 
phase gate is illustrated in Fig. El Each of the general- 
ized rotations that make up the gate is directly associated 
with of one of the measurements made in the eigenbasis 
of C/ 2 [±0/4]cr x C/ z [±0/4]L An initial cluster-state corre- 
lations which is used for the realization of a generalized 
rotation is shown in Fig. ^| the rotation uit 1 '^ [0/4] is 
realized via the measurement of the cluster qubit at the 
lattice site (3, 1) in the appropriate basis. 

The sign factors of the angles that specify the measure- 
ment bases depend on the outcome of a^-measurements 
only. Thus, after all a x - measurements have been per- 
formed, the measurement bases for the remaining qubits 
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can be deduced and the Toffoli phase-gate is realized in a 
single further time-step. The measurement pattern real- 
izes the generalized rotations directly and is not derived 
from a quantum logic network. 

Now we describe the realization of a four-qubit gate 
CARRY, which has one target and three control qubits. 
It performs a phase-flip a z on the target if at least two 
of the control qubits are in state |1) and otherwise does 



nothing, i.e. 



U CARRY = CXp -27T 



i=000 d I w(i)>2 



(175) 

Expanding the projectors on the control qubits into prod- 
ucts of Pauli operators one obtains 



J 



U, 



CARRY 



exp (-^«4 C3) ) exp (-i^crfcO) e xp (i^) exp exp 



exp 



(i J«r«) exp (-i^M^) exp (-i^fcO^W) exp (i^)^)) 



C/i, 



r 



(176) 
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FIG. 19: The three qubit controlled gate. Qubits displayed 
as squares in light gray are measured in the cr^-eigenbasis, the 
qubit displayed in dark gray is measured in the o-^-eigenbasis, 
and the measurement bases of the qubits displayed as framed 
squares are adaptive. 



The global phase is henceforth discarded. 

The unitary transformation is now subdivided into two 
parts, 



U CARRY = Uh,i U a -c 



(177) 



with U a - g = U g UfU e U d U c U b U a and U h4 = U l U h . Cor- 
respondingly, the cluster on which V carry is realized 
is divided into two sub-clusters. On the first sub-cluster 
the transformations U a to U g are realized, on the sec- 
ond sub-cluster Uh,i- The measurement pattern to real- 
ize Ucarry is displayed in Fie. 1191 The first sub-cluster 
stretches from x = to x = 8, with the input at x = 
and the intermediate output at x = 8. The qubits with 
8 < x < 16 belong to the second sub-cluster. 

Let us now explain the sub-gate U a - g - The conversion 
of the sequence <|176[) of generalized rotations is as in the 
previous examples. For each generalized rotation there 
is one cluster qubit in CM(U a - g ) whose measurement ba- 
sis specifies the respective rotation angle. Specifically, 
the measurement of the cluster qubit (3, 4) sets the ro- 



tation angle of U a , the measurement of qubit (4,3) sets 
the angle for {/&, (5,6) sets U c , (6,7) sets Ud, (6,5) sets 
U e , (6,3) sets Uf and qubit (6,1) sets U g . The quan- 
tum correlations of the initial cluster state which induce 
via the measurements of the cluster qubits in Cjf(i7 _ fl ) 
the quantum correlations associated with the generalized 
rotations are displayed in Fig. [3U1 

The realization of the gate requires two measurement 
rounds. In the first round the standard measurements of 
a x and a y are performed. Note that the rotation angle 
of Ud is twice as big as for the other rotations. To realize 
Ud of the cluster qubit (6, 7) the observable 



U z 



±- 



&xU z 



±er„ 



(178) 



is measured. Thus, the realization of Ud belongs to the 
first round of measurements. Strictly speaking, this mea- 
surement round does not belong to the gate but to the 
circuit as a whole since all standard measurements are 
performed simultaneously. 

In the second measurement round, of the remaining 
qubits in Cm(I7o-j) one measures the observables 



7T 
T77 



(179) 



The procedure to infer the sign factors in (|179fl and i|178|) 
is explained in Section III Fl 

The reason why the measurements in the tilted bases 
may all be performed simultaneously in the second round 
can be seen as follows. Be Qy the set of qubits measured 
in tilted bases. The contribution Uy,.q^ of the cluster 
qubits measured in tilted bases to the byproduct operator 
Ut, in (|84J) contains only a z-part but no a;-part. That is, 
it has the form 



(180) 



i£/C{t,Ci,C2,C3} 
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FIG. 20: Quantum correlations of the initial cluster state 
\tj>)c(u a - ) on the cluster C(U a - g )- These correlations induce 
via the o-^-measurements the quantum correlations for the 
state \ip') which act only on the output qubits and one cluster 
qubit in CM(U a - s )- The pattern of correlation centers in a) 
displays the correlation required to realize U a ; b),c) and d) 
display the correlations for Ub, U c and (7 e , respectively. The 
correlations used for the realization of Ud, Uf and U g are not 
shown. They are analogous to the one in d) used for the 
realization of U e - 
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FIG. 21: Quantum correlations of the initial cluster states 
on C(Uh) and C(Ui). These correlations induce, via the 
aaj-measurements, the quantum correlations for the states 
W)c(u h ) and \tl}')c(Ui) that involve only the respective out- 
put qubits and one qubit in the gate body. The pattern of 
correlation centers in a) displays the correlation required to 
realize Uh and b) the correlation for Ui. 





FIG. 22: In the three-qubit controlled gate CARRY, the tar- 
get qubit may travel either back or forth. 



In 1|83|) the byproduct operator appears "on the wrong 
side" of U a -g as does the contribution U^,Qy When 
the order of the gate and the byproduct operator is ex- 
changed, the byproduct operator may modify the gate. 
While this is, not surprisingly, indeed the case for the 
whole Us, it is not so for the contribution Uy,.q^ com- 
ing from the measurements in the tilted bases. Because 
Ut,.q^ bas only a z-part it commutes with U a — g . There- 
fore, the results of measurements in a tilted basis do not 
mutually affect the choice of their measurement bases. 

The fact that hat the byproduct operator U-^^q^, is 
indeed of form (|18Ufl we do not show here explicitly. For 
the byproduct operator created in the measurement of 
qubit (3, 4) realizing the transformation U a it may be 
verified from equation (|157fl in Section HV Bl 

The explanation of the second sub-gate, Uh,i, is anal- 
ogous. Fig. 1211 displays the quantum correlations of 
the initial cluster state which, via the measurements in 
CM(Uh,i), induce the required quantum correlations as- 
sociated with Ut and Ui. 

Two further points we would like to address in this sec- 
tion. The first is to note that the whole gate U carry can 
be performed on the QC C in two measurement rounds. 
The first measurement round is that of the a x -, a y - and 
a z -measurements which, strictly speaking, does not be- 
long to the gate but to the circuit as a whole. The second 



measurement round is that of the simultaneous measure- 
ments in tilted measurement bases. 

We have already seen that the measurements that real- 
ize the unitary transformations U a , ■■ , U g may be realized 
simultaneously, and this argument may be extended to 
the entire gate U carry- All the byproduct operators 
created with the measurements in tilted bases have only 
a z- but no x-part. Therefore they all commute with 
U carry- Thus, to choose the right measurement bases 
neither of the measurements in a tilted basis that realizes 
one of the rotations U a , .. , f/, needs to wait for another 
measurement in a tilted basis. 

Second, note that the for Ucarry the target-input 
and the target-output can be interchanged, see Fig. 
This holds because the (conditional) phase-flip on the 
target qubit is its own inverse. Thus, the target qubit 
may travel through the gate backwards. This property 
also holds for the Toffoli phase gate. We will make use of 
it in the construction of the quantum adder in the next 
section. 



G. Circuit for addition 

The QC c -version of the quantum adder corresponds to 
the quantum logic network as given in [2j|, see Fig. 1231 
In this paper we use the three-qubit controlled phase gate 
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FIG. 23: Quantum logic network for 4-qubit adder, c — a + 
b mod2 4 . The adder network is taken from I2H . The two- 
qubit controlled gate in this network is the Tofolli phase gate 
as discussed in Section ftV Fl A straightforward simulation of 
this network on the QC C would result in a quadratic scaling 
of spatial resources. However, the more compact realization 
discussed below requires only a linear overhead. 



CARRY together with a prior and subsequent Hadamard 
gate on the target qubit while in [23| the equivalent three- 
qubit controlled spin-flip gate is used directly. 

At first sight it appears as if the horizontal dimension 
of the cluster to realize the adder circuit would grow lin- 
early with the number of logical qubits n. This is, how- 
ever, not the case. The QC c -circuit may be formed in 
such a way that the horizontal size of the required clus- 
ter is constant such that the cluster size increases only 
linearly with the number n of logical qubits. To see what 
the QC c -realization of the quantum adder will look like, 
the network displayed in Fig. 1231 may be bent in a way 
displayed in Fig. |2U 

To "bend a network" is a rather informal notion. We 
therefore now specify what we mean by this. If a quan- 
tum circuit is displayed as a quantum logic network, the 
vertical axis usually denotes some spatial dimension, i.e. 
the location of the qubit carriers, and the horizontal axis 
corresponds to the sequence of steps of a quantum com- 
putation, i.e. a logical time. As the basic blocks of quan- 
tum computation in the network model, the universal 
gates, are unitary transformations generated by suitably 
chosen Hamiltonians the logical time becomes associated 
with physical time. This is, however, a peculiarity of 
the network model. If on the QC C a quantum logic net- 
work is simulated, the temporal axis is converted into 
an additional spatial axis. The temporal axis in a QC C - 
computation emerges anew. It has no counterpart in the 
network model. If we modify a quantum logic network 
in such a way that qubits travel from right to left, as 
done in Fig. [2] it does not mean that we propose to 
use particles that travel backwards in time because we 
do not need to respect the temporal axis implied by the 
network model. If one wants a semi-network picture that 
accounts for this, one may imagine the logical qubits as 
traveling through pipes on a two-dimensional surface. 

The reason why we may let the auxiliary qubits travel 
"backwards" is the identity displayed in Fig. [221 This 
arrangement of gates makes the circuit more compact. 
To complete the description of components from which 
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FIG. 24: Quantum logic network for 4-qubit adder, bent. 



a) 



b) 



■ 
■ ■ 



FIG. 25: Combination of two CNOT-gates (a) and its QCc- 
realization (b). 



the QC c -version of the quantum adder is built, a compact 
measurement pattern for the two combined CNOT-gates 
is displayed in Fig. 12*51 

The realization of the quantum adder in the network 
layout of Fig. |^| directly leads to the QC C -circuit for 
the quantum adder displayed in Fig. Please note 
that the displayed QC C -adder is for eight qubits while 
the networks in Figs. 1231 and 1241 are only for four qubits. 

For the quantum adder circuit in Fig. 1261 we have made 
two further minor simplifications. The first concerns the 
ancilla preparation. To prepare an ancilla qubit on the 
cluster in the state |+) means to measure the respective 
cluster qubit in the er^-cigenbasis (the randomness of the 
measurement outcome does not jeopardize the determin- 
istic character of the circuit). As can be seen from the 
Toffoli gate and the three-qubit controlled gate displayed 
in Figs. 1171 and !19l the ancilla qubits are located on clus- 
ter qubits which have only one next neighbor. As can be 
verified from the eigenvalue equations Q, to measure a 
qubit of a cluster state which only has one next neighbor 
in the eigenbasis of a x also has the effect of projecting 
this neighboring cluster qubit into an eigenstate of a z . 
Such cluster qubits may be removed from the cluster as 
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FIG. 26: Quantum adding circuit for two 8-qubit states. As in all figures displaying QC c -circuits, squares 
in light and dark gray represent cluster qubits measured in the a x - and <j y -eigenbasis, respectively. The 
measurement bases of qubits displayed as framed squares are adaptive. 



explained in Section Hi CI With these neighboring qubits 
removed the cluster qubits on which the initial ancilla 
qubits were located become disconnected from the re- 
maining cluster and may thus be removed as well. With 
the same argument, the cluster qubits carrying the an- 
cillas in their output state, and their next neighbors may 
also be removed. 

Second, between the QC c -realization of the CARRY- 
gates on the left and the subsequent blocks of CNOT- 
gates we have removed pairs of adjacent cluster qubits 
that would be measured in the eigenbasis of a x . Why 
this can be done has been explained for adjacent qubits 
in wires in Section III G 21 Here the situation is a little 
more involved since, like in case of the circuit for Fourier 
transformation displayed in Section llV El one of the re- 



moved qubits in each pair has more than two neighbors. 
But the method still works as can be easily verified. 

Let us now briefly discuss the resources required for 
the QC C -realization of an ?i-qubit adder. As can be seen 
directly from the circuit displayed in Fig. 1261 and the un- 
derlying network shown in Fig. 1241 with its repeating sub- 
structure, the adder requires a cluster of height 8n— 5 and 
of constant width 38. Thus the spatial and operational 
resources are, to leading order, 

S = = 304 n. (181) 

Concerning the temporal resources note that each pair 
of three-qubit controlled phase gates using the same con- 
trol qubits and the pair of Toffoli phase gates may be 
completed at one time instant but that one pair of gates 
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is completed after another. The reason why the mea- 
surements in the tilted bases that complete each pair of 
gates may be performed simultaneously is the same as 
the one given previously for the measurements in tilted 
bases of a single three-qubit controlled gate. The propa- 
gation of byproduct operators is most easily followed in 
the network of Fig. EH The temporal complexity T of an 
n-qubit QC c -addcr is 

T = n, (182) 

plus one step of a x -,a y - and cjj -measurements for the 
entire circuit. 

The corresponding network resources are to leading or- 
der S q \ n = 3n and O q i n = T q i n = 8n. For the count- 
ing of the operational and temporal network resources 
we have assumed that the three-qubit controlled spin- 
flip gate used in the addition circuit is composed of two 
Toffoli gates and one CNOT-gate as described in f23| . 
and that the CNOT- and the Toffoli-gate are regarded as 
elementary. 

Thus we find for both the network and the QC C - 
realization of the quantum adder that the spatial, tem- 
poral and operational resources scale linearly with n. 
Therefore, the resource overheads in one realization as 
compared to the other one are only constant. For the 
QC C this is much better than what is indicated by the 
bounds (|138[) . (|139[) and (|14L)|) . in particular for the spa- 
tial and operational resources. Equation (|139|l yields an 
upper bound on S which is ~ n 3 and (|140J) gives bounds 
on O and S which are ~ n 2 . Thus, the quantum adder 
is an example for which these bounds are very loose. In 
general they should not be mistaken as estimates. 

If the pre-factors are compared, one finds that for the 
realization of a quantum adder the QC C requires about 
100 times more spatial and 38 times more operational 
resources, while it is 8 times faster. However, since we 
compare different objects these ratios do not mean much 
apart from the fact that they are constant. It may be ar- 
gued that in case of the QC C spatial resources are not as 
precious as they usually are, for to create cluster states 
one needs a system with non-selective uniform interac- 
tion only while for quantum logic networks one gener- 
ally requires a system with selective interactions among 
the qubits. Concerning the operational and temporal 
resources, the QC C only uses one-qubit measurements 
while the corresponding network uses two- and three- 
qubit gates as elementary operations. 

H. Remarks 

We would like to add two remarks, one with regard to 
the elementary constituents of the QC C , and one with 
regard to their composition principle. 

For the particular set of gate simulations used in the 
QC C universality proof in Section^ the CNOT-gate and 
arbitrary one-qubit rotations, there is only a single in- 
stance of where one of these gates has been used as part 



of a more complicated gate in all examples of this sec- 
tion. Namely, the next-neighbor CNOT-gate has been 
used as part of the long-distance CNOT described in 
Section IIV CI Of universal gate simulations one might 
expect that any circuit is composed of them rather that 
they occur almost not at all. One could say, though, 
that the used set of gates is not a good choice for the 
universal set. In fact, in realizations of network quantum 
computers it is often the physics of the specific imple- 
mentation that determines which gates are elementary. 
For the QC C this is not so. The QC C may simulate, for 
example, general one-qubit rotations and Toffoli gates 
alike. Any gate simulation may be called "elementary" 
with the same right as any other, but they cannot be all 
elementary. The elementary constituents of the QC C are 
not gate simulations. 

As a consequence, the composition principle for these 
elements will be different from gate composition. At first 
sight, if we go through the examples of this section, we 
find that this is not yet reflected in the larger and more 
complicated constructions. For the quantum Fourier 
transform and the addition circuit we have, though play- 
ing with some tricks, ultimately imitated network com- 
position. 

However, in the smaller gates and sub-circuits such as 
the controlled phase gate, the Toffoli phase gate and the 
gate CARRY we find something that might give rise to a 
new and more appropriate composition principle. First, 
for the QC C it is not the one-qubit and two-qubit opera- 
tions that are particularly simple. In the Hamiltonian 
simulation circuit of Section IIV Bl wc found that it is 
easy to realize generalized rotations exp(iip a^ J ') where 

cr( J ) is a composite Pauli operator, = a6J (r[°', 
k a = x,y, z. Furthermore, in the subsequent examples 
of the multi-qubit gates in Sections IIVDI and IIV Fl we 
have decomposed the gates into such generalized rota- 
tions rather than into known standard gates on fewer 
qubits. 

Any unitary transformation may be decomposed into 
a unitary transformation in the Clifford group followed 
by generalized rotations. So, is this a new composition 
principle? With our present state of knowledge, the an- 
swer must be "Not yet." . First, though any transforma- 
tion may be rewritten in this form, it is presently not 
clear how to design quantum algorithms with these el- 
ements directly. Second, the construction uses the very 
concept of applying unitary transformations to the state 
of a quantum register. However, as we have explained in 
and also briefly sketched in Section llll Al the QC C has 
no quantum register. So, the generalized rotations and 
their concatenation at least have to be reformulated to fit 
the description of the QC C . In particular, they have to be 
made compatible with the graph states identified in Sec- 
tion llll Bi as characteristic quantum resource to represent 
algorithms. Nevertheless, it appears that the generalized 
rotations should be reflected in what may emerge as el- 
ementary constituents and composition principle for the 
QC C . 
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V. COMPUTATION WITH LIMITED SPATIAL 
RESOURCES AND IN THE PRESENCE OF 
DECOHERENCE 

In this section we describe how to perform QC C - 
computation on finite and possibly small clusters. If the 
cluster that may be provided by a specific device is too 
small for a certain measurement pattern it does not mean 
that the respective QC c -algorithm cannot be run on this 
device. Instead, the QC c -computation may be split into 
several parts such that each of those parts fits on the 
cluster. 

To see this consider Scheme ^ for the realization of 
gates. Scheme ^is applicable to any gate or sub-circuit. 
It is thus possible to divide the circuit into sub-circuits 
each of which fits onto the cluster. The adapted scheme 
is a process of repetitive re-entangling steps alternating 
with rounds of measurements. 

Specifically, one starts with the realization of the first 
sub-circuit acting on the fiducial input state located on 
h CC. The fiducial input is, while being processed, tele- 
ported to some subset 0\ of the cluster C. The set 0\ 
of qubits forms the intermediate output of the first sub- 
circuit. These qubits remain unmeasured while all the 
other qubits are measured to realize the first sub-circuit. 
Now the realization of the second sub-circuit begins. Its 
input state has already been prepared, 1% = 0\. The 
cluster qubits a € C\0± which have been measured in 
the realization of the first sub-circuit are now prepared 
individually in the state \+) a - This completes step 1 of 
Scheme n to realize the second sub-circuit. Step 2 is to 
entangle the whole cluster via the Ising interaction. In 
the third step all cluster qubits except those of the inter- 
mediate output O2 are measured whereby the realization 
of the second sub-circuit is completed. The intermediate 
output is now located on Oi. For the realization of the 
subsequent sub-circuits one proceeds accordingly. 

An advantage of this modified procedure is that one 
gets by with smaller clusters. A disadvantage is that the 
Clifford part of the circuit may no longer be performed 
in a single time step. 

Perhaps the most important advantage of the above 
construction is that in this way a basic requirement to 
make the QC C fault-tolerant can be fulfilled. Namely, 
decoherence can be controlled. If a single large clus- 
ter is used the computation might reach certain cluster 
qubits only after a long time such that the cluster would 
have already decohered significantly and it is not clear 
how error-correction could help in such a situation. This 
might, for any error rate, limit the duration of a com- 
putation. On the contrary, if the computation is split 
then the size of the sub-circuits may be adjusted such 
that each of them can be performed within a fixed time 
T and in this way, each cluster qubit is, before being 
measured, exposed to a bounded amount of decoherence 
specified by T. Thus, "fresh" qubits for computation are 
always provided. 



VI. CONCLUSION 

In this paper we have given a detailed account of the 
one-way quantum computer. We have shown that the 
QC C can be regarded as a simulator of quantum logic net- 
works. This way, we clarified the relation of the QC C to 
the network model of quantum computation and gave the 
universality proof. 

We have based our description on the correlations ex- 
hibited by cluster states, and states that can be created 
from them under one-qubit measurements. For this pur- 
pose, theorem 1 of Section III Fl is an important tool. It 
relates unitary transformations to quantum correlations 
exhibited by certain pure states. 

In Section llVI we have presented a number of example 
circuits such as the circuit for quantum Fourier transfor- 
mation and for addition. In this way, hopefully, we also 
have acquainted the reader with a number of construction 
principles for QC c -circuits. Note that the simulations of 
the universal gates required in the universality proof are 
hardly used. Instead, more compact measurement pat- 
terns have been found. 

The main purpose of this paper has been to provide a 
comprehensive description of the QC C from the network 
perspective. Beyond that, we have pointed out the non- 
network aspects of the QCg, such as the different nature 
of information processing [3, E| , and the connection to 
mathematical graphs. 
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APPENDIX A: CLUSTER DECOMPOSITION 

In Section III Al we have associated the cluster state 
\4>{k\) r on a cluster Cm with a graph G(Cn, Eq n ) where 
the set Ec N of edges is defined in the same way as in i|22|) 
for Ec ■ To decompose the cluster means in more precise 
terms to decompose the associated graph G(Cn , Ec N ) 
into subgraphs, that is we decompose both the set of 
vertices, Cn, and the set of edges, Eq n . As in the 
set of vertices is decomposed into the subsets C(gi), the 
sub-clusters corresponding to the gates gi, 

m 

C N = \JC( 9i ). 

i=l 

Herein, the sets C(gi) are overlapping. The output ver- 
tices of some sub-cluster C(gk), corresponding to the out- 
put cluster qubits of the gate gu , form if they are not the 
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output vertices of the whole graph- (some of) the input 
vertices of other sub-clusters {£(<?/)}. We define the sets 
/ and O of input and output vertices of Cat, and the set 
of overlapping vertices Vj/o as follows: 

/ = {a E C N : 3i\a E C/(&) A ~^3j | a E C {gj)} 
O = {a e C N : Si \ a g Cj(ffi) A 3 j \ a g C (gj)} 



I/O 



= {a e Cn ■ 3i | a € C/(.g 4 ) A 3j \ a e Co(ffj)}- 



(Al) 

In the same way as we decompose the set of vertices, Cn, 
into subsets C (<?;); according to (|49(l . we decompose the 
set of edges, Ec N , into subsets E(gi), 



|JV| 

\jE( 9i ). 



Now, for the decomposition to be useful, the subsets 
C(gi) and E(gi) must fulfill a number of constraints. The 
first of them is that each pair (C(gi), E(gi)) is again a 
graph, G{C{gi),E{gi)). As in (|51(l . this requires in par- 
ticular, that the endpoints of all the edges in E{gt) are 
in C(gi), 



a) 



b) 



Out 




Out 




C x (g) 



CJ9) C a (g) 



CJgJ Y,o ' 'C (gJ 
C M (9j C M lgJ 



FIG. 27: Two forbidden divisions of a cluster into sub- 
clusters. The functioning of the 4-qubit gate in a) cannot 
be understood from the functioning of the two 2-qubit gates 
on the sub-clusters. The displayed subdivision is excluded by 
constraint 15111 . The situation in b) can be treated with a 
careful assignment of the edges to the subgraphs. In this pa- 
per, however, it is excluded for simplicity by the assignment 
1A2II together with the constraint (I A3I) . 



Va g C N | (3e g E(g l ) s.th. a g e) : a G C(gi). 



we require that 



This already excludes a cluster Cn as displayed in 
Fig. 127b. There, the cluster Cn is decomposed into sub- 
clusters C(gi) and C(g2)- But there are edges, namely 
those which connect C{g2) and C(g\), which can nei- 
ther be included in G{C{gi),E{g 1 )) nor G(C(g 2 ),E(g 2 )). 
Therefore, condition l|49|) cannot be satisfied and conse- 
quently the decomposition of Cat into C(gi) and C{g2) is 
not allowed. 

It is necessary to exclude a decomposition as in Figl27b 
as the circuit displayed there cannot be understood from 
its components on the sub-clusters C{g\) and C(<?2)- The 
reason for this is that the cluster states on the sub- 
clusters C{gi) and £"(.92) are mutually entangled. This 
is caused by precisely those interactions S ab that corre- 
spond to the edges in Ec N which could not be included 
in either of the subgraphs. 

The central condition (|5U|I is that, in contrast to the 
sets of vertices C(gi), the sets of edges E(gi) do not over- 
lap, 

Vi,j = l..lAT\,i?j: %)n%) = 0. 

We need a way to assign the edges to the vertices in 
the sub-clusters. In this paper, for simplicity we adopt 
the convention that sets of edges E{gi) are chosen such 
that the subgraphs G(C(gi), E{gij) are induced subgraphs 
in G(C n ,E Cn ), i.e. 

G(C{g i ),E{g i ))=G[C(g i )\. (A2) 

As we have overlapping vertices, this simple assignment 
may run into conflict with condition l|5(J|) . To avoid this, 



G[Ci( 9i )] = G(C/( ffi )),0), and 

G[C ( 9l )} = G(C (<?,)), 0), Vi = l..\Af\. {A6 > 

That is, the vertices within the sets Ci(gi), Co(gi) are 
not connected by edges. 

The assignment convention (|A2|) together with the con- 
dition (| A3|> exclude decompositions of a cluster like the 
one displayed in Fig 127b . In principle, however, such a 
decomposition is possible. In such a case the assignment 
of the edges to the sub-clusters is not as simple as IA2I . 
Instead, each edge had to be assigned to only one sub- 
cluster while respecting condition 151|) . In order to keep 
the notation as simple as possible, in this paper we do not 
consider decompositions of the sort displayed in Figl27b. 

The constraints displayed so far applied to the sub- 
graphs. We have two further constraints on the graph 
G(Cn, Eq n ), with regard to its decomposition. Or, to 
put it in different terms, the remaining constraints are 
for the proper composition of the graph out of the sub- 
graphs. They correspond to the usual gate composition 
rules. We require that each input and output vertex is 
the input or output vertex of exactly one gate, 

Va g /UVf/o : -ai,j\i^j, ae Cj{gi) A a € 
Va g O U Vf/ : ~^3i, j | i ^ j, a g C {gi) A a g C (g 3 )- 

(A4) 

Further, the graph G(Cn , Ec N ) with the vertices of / U 
OUVj/o and associated edges removed disintegrates into 
mutually disconnected induced subgraphs corresponding 
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to the gate bodies, 

G(C Ni E Cn )\(IU0UV i/o ) = \jG[C M ( 9i )]. (A5) 

i=l 

To summarize, the two central conditions and l|5Tj|) 
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